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ABSTRACT 

We measure and analyse the clustering of the Baryon Oscillation Spectroscopic 
Survey (BOSS) relative to the line-of-sight (LOS), for LOWZ and CMASS galaxy 
samples drawn from the final Data Release 12 (DR12). The LOWZ sample contains 
361 762 galaxies with an effective redshift of ziowz = 0.32, and the CMASS sam¬ 
ple 777202 galaxies with an effective redshift of ^cmass = 0.57. From the power 
spectrum monopole and quadrupole moments aronnd the LOS, we measure the 
growth of structure parameter / times the amplitude of dark matter density flnc- 
tuations as by modeling the Redshift-Space Distortion signal. When the geometrical 
Alcock-Paczynski effect is also constrained from the same data, we find joint con¬ 
straints on /as, the prodnct of the Hubble constant and the comoving sound hori¬ 
zon at the baryon drag epoch H{z)rs{zd)^ and the angular distance parameter di¬ 
vided by the sound horizon DA{z)lrs{zd)- We find /(-2iowz)^s('2^iowz) = 0.394 zb 0.062, 
DA{z\ow7)lrs{zd) = 6.35 ± 0.19, H{ziowz)rsizd) = (11.41 ± 0.56) 10^kms“i for the 
LOWZ sample, and /(zcmass)cr8(2cmass) = 0.444 ± 0.038, L>yi(zcmass)/rs(zd) = 9.42 ± 
0.15, H{zcmass)rsizd) = (13.92 ± 0.44) lO^kms”^ for the CMASS sample. We find 
general agreement with previous BOSS DRll measurements. Assuming the Hubble 
parameter and angular distance parameter are fixed at fiducial ACDM values, we 
find /(ziowz)o-8(Aowz) = 0.485 ± 0.044 and f{zcniass)(^8izcmass) = 0.436 ± 0.022 for the 
LOWZ and CMASS samples, respectively. 

Key words: cosmology: theory - cosmology: cosmological parameters - cosmology: 
large-scale structure of Universe - galaxies: haloes 
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1 INTRODUCTION 

The large-scale distribution of matter, as observed through 
galaxy clustering, encodes significant cosmological informa¬ 
tion. Much of this can be extracted from the shape and am¬ 
plitude of the galaxy power spectrum multipole moments. 
The focus of this paper is the observed anisotropic redshift- 
space distortions (RSD; [Kaiser||19^ I caused by peculiar ve¬ 
locities, which contain information about how gravity be¬ 
haves at large scales and about the total matter content 
of the Universe. As these distortions depend on the growth 
of structure, they offer an independent and complementary 
technique to measure the matter content and to test gravity, 
compared to studies of the cosmic expansion history. 

In this paper we measure the galaxy power spectrum 
monopole and quadrupole moments calculated from the 


galaxy sample of the Sloan Digital Sky Survey III ( Eisen- 
stein et al. 201 1| Baryon Oscillation Spectroscopic Survey 


(BOSS; [Bolton et al.||20l^ [Dawson et al.||2013| jSmee et al.j 

2013) data release 12 (DR12; Alam et al. 2015a I. The galaxy 
catalogues drawn from the final data release from BOSS, 
DR12, cover the largest cosmic volu me ever observed , with 
an effective volume 14ff = 7.4Gpc^ ( Reid et al.||2015 |. The 
number of independent modes contained allows us to ob¬ 
serve the RSD with the highest ever significance, and use 
them to set tight constraints on the growth of structure / 
times the amplitude of primordial dark matter power spec¬ 
trum erg at the effective redshifts of 2 = 0.32 and 2 = 0.57. 
Additionally we are able to constrain DA{zes)/rs{zd) and 
H{zefi)rs{zd) by performing the Alcock-Paezynski test (AP; 
Alcock fc Paczynski||1979 ). 

This paper forms part of an initial set of DR12 papers 
produced by the BOSS galaxy clustering team. In ( jCuestaj 
et al. 2015 companion paper) we measure BAO from the 


correlation function moments measuring the cosmological 
expansion history. A complimentary approach is provided 
in ( Gil-Marm et al.]|2015b companion paper), which looks 
at the Baryon Acoustic Oscillation (BAO) in power spec¬ 
trum monopole and quadrupole moments, with both papers 
applying an algorithm to reconstruct the primordial density 
distribution. The BOSS galaxy targeting and catalogue cre¬ 
ation algorithms are presented in [Reid et al.j ( |2015| ), which 
includes an extensive analysis of the catalogues themselves, 
and the methods employed to correct for observational ef¬ 
fects. 

RSD measurements have a long history following the 


seminal theoretical paper of Kaiser (19871, and the develop¬ 


ment of large galaxy surveys. Key milestones include mea¬ 


surements from the PSCz galaxy survey (Tadros et al. 19991, 


2-degree Field Galaxy Redshift Survey (Peacock et al.|2001 


Hawkins et al. 2003 Percival et al. 20041, WiggleZ ([Blake 


et al. 2012), 6-degree Field Galaxy Survey (Beutler et al. 


2012) and the SDSS-II LRGs (Oka et al. 2014). Previous 


measurements made from BOSS including DR9, when the 


survey was approximately one third complete (Samushia 
e£^alj20^|^ and most recently from DRlO/11 IChuang^e^^ 


2013 Beutler et al.|2014l[Samushia et al.|2014|[Sanchez et al. 

2014 Reid et al.[[2014 Alam et al.pOlhb ). 

The RSD measurements made from the BOSS DRll 
sample, which are discussed further in 0 exhibit a scatter 
that is not negligible with respect to the statistical errors. 
This is caused by adopting different approaches to model 


quasi-linear and non-linear behaviour in the measured clus¬ 
tering, beyond the linear RSD signal. As the signal-to-noise 
increases as we move to small scales, RSD measurements are 
very sensitive to the behaviour of the model on the small¬ 
est scales fitted. It is therefore important that the model 
is general enough to encompass all of the unknown be¬ 
haviour, whilst also being able to accurately model the data 
itself. Gonsequently, a large effort on developing a model 
that is able to account for: i) mode-coupling; ii) galaxy 
bias; iii) RSD ; up to quasi- and non-linear scales is re¬ 
quired. Some improvements have recently been achieved on 
i) beyond the linear regime for the real space power spec¬ 
trum using perturbation theory schemes. There are several 
models that attempt to do this task: cosmological standard 


perturbation theory (see Bernardeau et al. 2002 and ref- 


erences therein), Lagrangian perturbation theories (Hivon 

et al.||1995 

|Matsubara||2008 [Gar Ison, White & Padm 

anab- 

ageas 

han 2009| 

Okamura, Taruya & Matsubara||2011| Va 

& Nishimichi 2011), time renormalization (Pietroni 

2008 

Anselmi & Pietroni 2012), Eulerian resumed perturbation 

theories (( 

I!rocce & Scoccimarro 2006 Bernardeau, Crocce 

& Scoccimarro||2008||Wang & Szalay|2012||Elia et al. 

20111 

Taruya et al. 2012 Bernardeau, Crocce & Scoccimarro 

20121 

and closure theory (Taruya & Hiramatsu|2008 ). Recent im- 


provements in the galaxy bias model describe accurately 
how the galaxies trace dark matter, including non-linear 


and non-local terms (Nishimichi & Taruya 2011 McDon- 
[ald k, Roy[[2009{ [Saito et al.[[2014|~ in addition to primor¬ 


dial non-Gaussian terms (Biagetti et al. 2014). The final 


point iii) is to accurately model the mapping from real to 
redshift space statistics. Different approaches to this prob¬ 


lem includes: the TNS model (Taruya, Nishimichi & Saito 


2012 


2010[), the Distribution function approach model (Vlah et al. 
[Okumura, Seljak fc Desjacques[|2012[) and the Gaus¬ 


sian streaming model (Reid & White 2011). 


In this paper we use a resumed perturbation theory 
approach (RPT) in order to describe the dark matter non¬ 
linear power spectra components as it is described in 


Marfn et al. (2012). We combine this approach with the non¬ 


linear and non-local galaxy bias model of[McDonald fc Roy 


(2009), as it was done in previous analyses (Beutler et al. 


2014 


Gil-Marm et al. 2015a). Finally, we account for the 


RSD using the redshift space mapping presente d in[Taruya,[ 
Nishimichi & Saito ( 2010| ; Nishimichi & Taruya (2011). This 


modelling has demonstrated to describe accurately the RSD 


for both dark matter particles and dark matter haloes (Gil- 
Marm et al. [2014|. Note that in previous approaches. 


Beutler et al. (2014) use a different model of RPT to de¬ 


scribe the dark matter components and a different fc-range 
is considered for obtaining the best-fit parameters. Also in 
this paper we treat / and erg as a free independent parame¬ 
ters during the fitting process, and is only at the end of the 
analysis which we constrain /erg. On the other hand [Beut-[ 


ler et al. (2014) constrain /erg fixing the value of erg on the 


non-linear components of the model, which do not depend 
on the particular combination / x erg. 

This paper is organised as follows. In !|^we present the 
description of the LOWZ and GMASS DR12 data samples 
and the resources used for computing the covariance matri¬ 
ces and for testing the theoretical models. In )|^we present 
the estimator used for measuring the power spectrum multi¬ 
poles. In ^the results including the best-fit parameters and 
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their errors are presented. i|^contains information about the 
theoretical models used to describe the galaxy power spec¬ 
trum multipoles in redshift space. In this section, we also 
include information about how to model the AP distortions 
and the effect of the survey window in the measurements. 
In we present the details about the parameter estima¬ 
tion, including how the covariance matrices are extracted 
and how the best-fit parameters and their errors have been 
computed. In we present a final analysis of the results 
presented in Cpmd how their compare with other galaxy 
surveys and other BOSS analyses, as well as the effect that 
changing the cosmological model has on /erg. Finally, in !|^ 
we present the conclusions of this paper. 


2 DATA AND MOCKS 


2.1 The SDSS III BOSS data 


As part of the Sloan Digital Sky Survey III (Eisenstein 
et al. I IMTT]) the Baryon Oscillations Spectroscopic Survey 


(BOSS) (Dawson et al. 20131 measured spectroscopic red- 
shifts for more than 1 million galaxies and over 200 000 
quasars. The galaxies were selected from multi-colour SDSS 
imaging (Fukugita et al.||1996 Gunn et al.|1998 Smith et al. 


|2002[ |Gunn et al.|2006[|Doi et al.|2010p focussing on the red- 
shift range of 0.15 SS z Si 0.70. The galaxy survey used two 
primary target algorithms, selecting samples called LOWZ, 
with 361 762 galaxies in the final data release DR12 (|Alam| 


et al. 2015a I between 0.15 ^ a si 0.43 and GMASS, with 


777 202 galaxies in DR12 between 0.43 ^ 2 ^ 0.70. The full 
targeting algorithms used and the method for calculating the 
galaxy and random catalogues are presented in |Reid et al.| 
(20151, which also shows that the samples jointly cover a 


large cosmic volume I4ff = 7.4 Gpe'^ with a number density 
of galaxies as a function of observed redshift, that ensures 
that the shot noise does not dominate at BAG scales. Full 
details of the catalogues are provided in |Reid et al.| ( |2015[ ), 
and we do not replicate this here. 

In order to correct for several observational artifacts 
in the catalogues, the GMASS and LOWZ samples incor¬ 
porate weights: a redshift failure weight, w^i, a fibre col¬ 
lision weight, Wfc, and a systematic weight, Wsys (GMASS 
only), which combines a seeing condition weight and a stel¬ 
lar weight (Ross et al.|2012 Anderson et al.|2014[|Reid et al. 


2015). Hence, each galaxy target contributes to our estimate 


of the target galaxy density field by 

Wq — -U^sys (rUrf -(- rufe 1). 


( 1 ) 


The redshift failure weights account for galaxies that have 
been observed, but whose redshifts have not been measured: 
nearby galaxies, which are approximated as being “equiva¬ 
lent” are up-weighted to remove any bias in the resulting 
held. The hbre collision weight similarly corrects for galax¬ 
ies that could not be observed as there was another tar¬ 
get within 62", a physical limitation of the spectrograph 
(see [Ross et al. (2012l for details). The systematic weight 
accounts for fluctuations in the target density caused by 
changes in the observational efficiency. It is only present for 
the GMASS sample, which relies on deeper imaging data, 
and such a weight is not required for the brighter LOWZ 
sample. 
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Additionally, we use the standard weight to balance re¬ 


gions of high and low density (Feldman, Kaiser & Peacock 
1994 Rentier et aI.|[2M4 l, 


WFKp(r) = 


.(r) 


Wsys(r) -f Wc(r)n(r)Pbao ’ 


( 2 ) 


where n is the mean number density of galaxies and Pbao 
is the amplitude of the galaxy power spectrum at the 
scale where the error is minimised. We assume Pbao = 
10 000 Mpc^/i“®, which corresponds to the ampl itude of the 
power spectrum at scales k ~ O.lOhMpc”^ (Reid et al. 


20151. 


2.2 The mock survey catalogues 

Galaxy mock catalogues have become an essential tool in the 
analysis of precision cosmological data provided by galaxy 
surveys. They provide a fundamental test of large-scale 
structure analyses and help to determine errors on measure¬ 
ments. As much of the large-scale physics can be captured 
using approximate methods, we do not necessarily need to 
base mock catalogues on full N-body cosmological simula¬ 
tions: structure formation models can be calibrated with 
a small number of N-body simulations, and the parameter 
space studied using a more efficient scheme. In this paper we 
use mocks cerated by two different approaches: MultiDark- 
Patchy BOSS DR12 mock^ (hereafter MD-Patchy) (Ki- 


taura et al. | [20T^ companion paper) and Quick-Particle- 


Mesh mocks (hereafter QPm) ([White, Tinker fc McBride 


20141. Both schemes incorporate observational effects in¬ 


cluding the survey selection window, veto mask and fiber 
collisions. 

MD-Patchy mocks rely on Augmented Lagrangian 
Perturbation Theory (ALPT) formalism ( jKitaura &: Hefij 
20131, which is based on splitting the displacement field 


into a long and a short-range component. The long-range 
component is computed by second order Lagrangian Pertur¬ 
bation Theory (2LPT), whereas the short-range component 
is modeled using the spherical collapse approximation. 
The MD-Patchy mocks use 10 combined snapshots at 
2 = 0.1885, 0.2702, 0.3153, 0.3581, 0.3922, 0.4656, 0.5053, 
0.5328, 0.5763, 0.6383. The underling cosmology for these 
mocks has been chosen to be (Ha, ^b, os, h, Us) = 
(0.692885, 0.307115, 0.048, 0.8288, 0.6777, 0.96), being 
very close to the best-ht values of the last release of 
PlancklB ( [Planck Collaboration et al.|[2015p . For this 
cosmological model the sound horizon at the baryon-drag 
redshift is rs{zd) = 147.66Mpc. 

The QPM mocks are based on low-resolution par¬ 
ticle mesh simulations that accurately reproduce the 
large-scale dark matter density field, in combination 
with the halo occupation distribution technique (HOD) 
to populate the resolved haloes with galaxies. For 
the QPM mocks, the snapshots are at the effective 
redshift of, 2 eff = 0.55 for GMASS and 2efi = 

0.40 for LOWZ. The underling cosmology for these 
mocks has been chosen to be (Ha, H^, H;,, as, h, Us) — 
(0.71, 0.29, 0.0458, 0.80, 0.7, 0.97). For this cosmological 


^ http://data.sdss3.org/datamodel/index-files.html 
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model the sound horizon at the baryon-drag redshift is 
raGd) = 147.10 Mpc. 

2.3 Fiducial Cosmology 

We have opted to analyse both mocks and data using the 
same cosmological model. The fiducial value assumed for 
this is = 0.31, which is in agreement with the last 
PlancklS release. As a consequence we will analyse the 
mocks using a value of that is different than their true 
values. When converting redshift into comoving distances, 
this will introduce an extra anisotropy to the one gener¬ 
ated by the peculiar velocities. In our analysis this effect 
will be accounted by the AP scaling relations presented 
in §5.3| The rest of cosmological parameters in the fidu¬ 
cial cosmology are n“) = 

(0.69, 0.31, 0.049, 0.8475, 0.6711, 0.9624). For this cosmol¬ 
ogy the sound horizon at the baryon-drag redshift is Vg (zd) = 
148.11 Mpc. 


3 MEASURING POWER SPECTRUM 
MOMENTS 


In order to compute the galaxy power spectrum we start 
by defining the Feldman-Kaiser-Peacock function ( [Feldman^ 
Kaiser fc Peacock|1994| ), 

K(r)w(r) - cm^r)], (3) 

■'2 

where n and Ug are, respectively, the observed number den¬ 
sity of galaxies and the number density of a synthetic catalog 
Poisson sampled with the same mask and selection func¬ 
tion as the survey with no other cosmological correlations. 
The functions Wc and wfkp were defined in Eqs and 
respectively. The factor a is the ratio between the weighted 
number of observed galaxies over the random catalogue, 
a = Wc/Ng, where Ng denotes the number of objects 

in the synthetic catalog and Agai the number of galaxies in 
the real catalog. The factor I 2 normalises the amplitude of 
the observed power in accordance with its definition in a 
universe with no survey selection. 


h = J <frwlKp{nWg)^{r). 


( 4 ) 


Following the Yamamoto estimator | Yamamoto et al.|2006 1, 
we define the multipole power spectrum estimator as. 




Sk) = 


{ 21 + 1 ) 


[I 

J dr2 F(r2)e"‘‘‘ ''"£r(k ■ fa) - 


.(k) 


( 5 ) 


where dflfe is the solid angle element, Ci is the Legendre 
polynomial of order RpoLson ^ke Poisson shot noise term, 

-PpoLon(k) = (1 + a) J dr h(r)w^(r)£r(k • f), (6) 

where the integral has been performed as a sum over the 
galaxy catalogue. For multipoles of order £ > 0, PpJ^^son ^ 
P^^\ and consequently the shot noise correction is negligi¬ 
ble. This estimator keeps the relevant LOS information by 


approximating the LOS of each pair of galaxies with the 
LOS of one of the two galaxies. This is a reliable approxi¬ 
mation on the scales of interest, which clearly improves on 
assuming a single fixed LOS for the whole survey for I > 0, 
but will eventually break down at very large scales and high 
order multipoles (Yoo & Seljak||2015 Samushia, Branchini 
fc Percival|M5| . 


The implementation of the Yamamoto estimator is per¬ 
formed using multiple FFTs (Fast Fourier Transform), each 
measuring the LOS-weighted clustering along different axes 
as presented in |Bianchi et al.| ( |2015[ ). Thus, the computation 
of the monopole and quadrupole can be written in terms of 
N log N processes (where N is the number of grid-cells used 
to discretise the galaxy field), which is significantly faster 
than performing the sum over galaxies used in previous anal¬ 
yses ( Beutler et al.|[2C)14 |. 

We use a random catalogue of number density of 
hs(r) = a“^h(r) with oT^ ~ 50. We place the LOWZ 
and CMASS galaxy samples on 1024® grids, of box side 
Li, = 2300/i“®Mpc for the LOWZ galaxies, and Li, = 
3500/i“®Mpc to fit the CMASS galaxies. This corresponds 
to a grid-cell resolution of 3.42/i“®Mpc for the CMASS 
galaxies and 2.25/i“®Mpc for the LOWZ galaxies. The 
fundamental wave-lengths are therefore kf = 1.795 • 

10“® /iMpc“® and kf — 2.732-10“® /iMpc“® for the CMASS 
and LOWZ galaxies, respectively. We have checked that for 
k ^ 0.3/iMpc“®, doubling the number of grid-cells per side, 
from 1024 to 2048, produces a negligible change in the power 
spectrum, <C 1%. This result indicates that using 1024® 
grid-cells provides sufficient resolution at the scales of in¬ 
terest. We apply the Cloud-in-Cells scheme (CiC) to as¬ 
sociate galaxies to grid-cells, and bin the power spectrum 
fc—modes in 60 bins between the fundamental frequency kf 
and a maximum frequency of ku = 0.33 LMpc, with width 
A login k = [logio(fcM) - logio(fc/)] /60. 


We limit the scales fitted as follows: our procedure 
for determining the largest scale we use for the fitting 
process is based on limiting the impact of the systematic 
weights, and is presented in Appendix]^ We limit scales to 
k > 0.02hMpc“® for the monopole and k > 0.04hMpc“® 
for the quadrupole. The smaller (larger) the minimum scale 
(fc-value) included, the more fc-modes are used and therefore 
the smaller the statistical errors of the estimated parame¬ 
ters. However, small scales are poorly modeled in compari¬ 
son to large scales, such that we expect the systematic errors 
to grow as the minimum scale decreases. Therefore, we em¬ 
pirically find a compromise between these two effects such 
that the systematic offset induced by poorly modelled non¬ 
linear behaviour is smaller than the statistical error. To do 
so, we perform different best-fit analysis for different mini¬ 
mum scales and check that the best-fit parameters of inter¬ 
est does not change significantly (compared to the statistical 
errors) as a function of this minimum scale. 


4 THE POWER SPECTRUM MULTIPOLES 

The top sub-panel of Fig. presents the power spectrum 
monopole (blue squares) and quadrupole (red circles) for 
LOWZ and CMASS DR12 data measurements (top and bot¬ 
tom panels as labeled) from the combination of the NGC 
and SGC data. This combination has performed by aver- 
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aging the NGC and SGC power spectra weighting by their 
effective area, 

pW _ (^P^f^Ai<iGC + -Pj§c^SGc)/(^NGC + ^SGC), (7) 

where Angc and Asgg are the effective areas of the NGG 
and SGG, respectively, whose values are ^ngc = 5836 deg^, 
= 2501 deg^ = 6851 deg" and = 

2525 deg^ The best-fit model predictions are shown by the 
solid lines, taking the average of the models fitted using 
the covariance extracted from QPM and MD-Patchy mocks 
with parameters as reported in Table in 0 Details about 
the models are presented in 0 and covariance matrices in 
§6.1[ The blue solid line shows the model for the monopole 
and the red solid line for the quadrupole. Error-bars cor¬ 
respond to the rms of the mocks, averaging between those 
calculated using covariance matrices determined with QPM 
or MD-Patchy mocks. The comparison of the rms of the 
QPM and MD-Patchy mocks can be found in Fig. and is 
discussed in § 6 . 1 | 

In the lower sub-panels, we present the fractional dif¬ 
ferences between the data and the best-fit model. For both 
LOWZ and CMASS samples, the model is able to repro¬ 
duce the monopole data points up to fc ~ 0.24/iMpc“^, 
within 3% accuracy (indicated by the black dotted horizon¬ 
tal lines). The model reproduces the measured quadrupole 
with an accuracy of ~ 10% for the LOWZ sample and of 
~ 5% for the GMASS sample. However, the signal-to-noise 
ratio of the observed quadrupole is not sufficiently high to 
determine if the observed fluctuations are statistical or sys- 
tematics of the model. For LOWZ and CMASS samples, the 
fitting process ignores the large scales {k ^ 0.02hMpc“^ 
for the monopole and k ^ 0.04/iMpc”'^ for the quadrupole) 
because of the effects of star contamination as it has been 
discussed in (2.1 (see Appendix [A] for a further discussion 
on how these limits have been decided). The effects of the 
fiber collisions on the best-fit parameters of the models are 
discussed in Appendix [B] 


5 MODELLING THE POWER SPECTRUM 
MULTIPOLES 

In this section we present the model used to analyse the 
monopole and quadrupole power spectra in 0 The model- 
lisation is done in the following four steps. 


(i) In 1 5.1 we present the galaxy bias model that maps the 


dark matter theoretical predictions into galaxy statistical 
observables. 


(ii) In [ 5.2 we present the model that relates the real 
space statistical moments with redshift-space ones. 


(iii) In 15.3 we incorporate AP parameters in order to 
allow changes due to inaccuracies when converting redshifts 
into distances by assuming a different value of than the 
actual one. 


^ The areas for both samples were initially planned to be the 
equal, as both samples are coming from the same spectroscopic 
pointing. However, due to difficulties during the early phases of 
the project, the sky area of the LOWZ sample lags that of t he 
full survey by approximately 1000deg" (Anderson et al.|2014l. 


(iv) In 15.4 we describe how the window survey mask is 
applied in order to account for observational effects due to 
the geometry of the survey. 


5.1 The bias model 

We assume an Eulerian non-linear and non-local bias model 


proposed by McDonald & Roy (20091 and previously used 


for analysing the power spectrum multipoles and bispec 
trum of DRll CMASS BOSS galaxies ( Beutler et al.|[^014 


Gil-Marfn et al.||2015a l. A priori the non-local galaxy bias 
model depends on 4 free parameters: the linear bias 6 i, the 
non-linear bias 62 and non-local bias parameters 6^2 and 
bsni- In order to reduce the number of free bias parame¬ 
ters, we assume that the bias is local in Lagrangian space, 
which sets the values of and fesni given the linear bias 
coefficient, bi. This assumption has been validated using N- 
body simulations, and it provides consistent results between 
the power spectrum and bispectrum for the CMASS sam¬ 
ple ( Gil-Marfn et al.]|2014 2015aI. In case the condition of 
local Lagrangian bias were relaxed, the parameters b ^2 and 
bani would be treated as free parameters, increasing, conse¬ 
quently, the number of free parameters of the bias model. 
In such case, the model would account for a more general 
galaxy biasing, but the error-bars on the parameters of in¬ 
terest would also increase. Since, the results of the power 
spectrum and bispectrum of mocks and N-body simulations 
suggest that the bias model of galaxies and haloes is consis¬ 
tent with local Lagrangian, we also assume it for the data. 


5.2 Modelling the redshift space distortions 

The mapping from real space to redshift space quantities 
involves the power spectrum of the velocity divergence. We 
assume that there is no velocity bias between the under¬ 
ling dark matter field and the galaxy field at least on the 
relatively large scales of interest. We follow the same red¬ 


shift space modelling that in previous analysis (Beutler 


et al. 2014 Gil-Marfn et al. 2015aI, described in Taruya, 


Nishimichi & Saito (20101 and Nishimichi & Taruya (20111, 


which provides a prediction for the redshift space power 
spectrum in terms of the matter-matter, velocity-velocity 
and matter-velocity non-linear power spectra. Expressions 
for the non-linear power spectra used in this paper, were 


presented in Gil-Marfn et al. (2015a I; the model for the non¬ 


linear matter quantities was obtained using resumed pertur¬ 
bation theory at 2-loop level as described in |Gil-Marm et ah] 
(20121. The necessary linear power spectrum input was com¬ 
puted using CAMB ( Lewis, Ghallinor fc Lasenby|200^ . 

We account for the Fingers-of-God (hereafter FoG), 
through the Lorentzian damping factor, as described in |Gil-| 
Marfn et al. (2015a I. This factor has one free parameter 


refered as upoG . With this parameter we aim to describe the 
non-linear damping due to the velocity dispersion of satellite 
galaxies inside host haloes. However we treat this factor as 
an effective parameter that encode our poor understanding 
of the non-linear RSD and we marginalised over. 

In this paper we consider that the shot noise contri¬ 
bution in the power spectrum monopole may be modified 
from that of a pure Poisson sampling. We parametrise this 
deviation through a free parameter, Anoise, i.e., = 
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Figure 1. The top sub-panels display the measured LOWZ- (top panel) and CMASS-DR12 (bottom panel), monopole (blue squares) 
and quadrupole (red circles) power spectra. For both cases, the measurements correspond to a combination of the northern and southern 
galaxy caps according to their effective areas as presented in Eqig The error-bars are the average values of the dispersion among 
realizations of the QPM and MD-Patchy mocks. The red and blue solid lines correspond to the best-fit model using the parameters listed 
in Table ^ for /cmax = 0.24/iMpc“^. For simplicity, we plot the average between models corresponding to the parameters obtained using 
QPM and MD-Patchy covariance matrices. The bottom sub-panels show the ratio between the power spectrum multipoles measurements 
and the best-fit model presented in the top sub-panel. The quadrupole symbols have been displaced horizontally for clarity. The dotted 
black lines represent a 3% deviation. 
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(1 - ^noise)Pp°Laon> w^ere the terms is the Pois¬ 

son predictions for the shot noise as is presented in Eq. 
For Anoise = 0 we recover the Poisson prediction, whereas 
when Anoise > 0 we obtain a sub-Poisson shot noise term 
and Anoise < 0 a super-Poisson noise term. 


5.3 The Alcock-Paczynski effect 


The AP effect (Alcock & Paczynski 19791 is caused by con¬ 


verting redshift into distance using a different cosmology 
from the actual one, which introduces a spurious anisotropy 
in the power spectrum that can be measured. Along the 
LOS, the observed signal is sensitive to the Hubble param¬ 
eter through oc H~^{z), when the clustering is measured on 
scales that are small compared with cosmological changes 
in the distance-redshift relationship. On the other hand, in 
the angular direction the distortions depend on the angular 
distance parameter, Da{z). When a fiducial model is as¬ 
sumed to convert redshifts into distances, the AP effect can 
be described by the dilation scales. 


a±_ = 


H{z)rs{zd) 

DA{z)Os^{Zd) 

Df{z)rs{zd) ’ 


( 8 ) 

(9) 


where and a± are the parallel- and perpendicular-to-the- 
LOS dilation scales, respectively. Here, and D^{z) 

are the fiducial values (those corresponding to the assumed 
cosmology to convert redshifts into distances) of the Hubble 
constant and the angular dimeter distance at a given redshift 
z, respectively. On the other hand, the fiducial sound horizon 
at the baryon-drag redshift is given by r^'^{zd)- The factors 
ay and a±_ describe how the true wave-length modes, fc'i and 
fcp, have been distorted into the observed ones, fey and fcx: 
fc|| = s-nd — a^kf, by the effect of assuming a 

different cosmological model. 

This component of observed anisotropy is modelled 
through the ay and a± parameters (see eq. 60 of 


Beut- 


ler et al. 20141. The geometric AP effect also affects the 


BAO scale, and the assumed distance-redshift relation has 
the potential to shift the BAO peak position differently in 
the monopole and quadrupole moments of the observed co¬ 
moving clustering signal. Therefore, the AP parameters are 
simultaneously measured in our analysis, even though our 
focus is on measuring the RSD. Analyses that do not wish 
to measure the RSD signal and focus on the BAO scale 
can make use of reconstruction and therefore, in general, 
provide better measurements of the AP effect (e.g. |Cuesta| 
et al.|20l5 Gil-Marm et al.|2015b I. We will comment on the 
potential differences further in Izl 

Assuming the fiducial cosmology described in §2.3[ 


the fiducial values for H{z) and Da{z) are, H" (^lo-wz) = 

999.23 Mpc for the 

and n (^cniass) “ 


DTizio^L = 

^lowz — 0.32, 

-^fid / 


79.49 kms“ Mpc” 

LOWZ sample at 

92.25kms“^Mpc“^, ^“(zcmass) = 1398.43Mpc for the 
CMASS sample at Zcmass = 0.57. The value for the fidu¬ 
cial sound horizon distance is r^‘^{zd) = 148.11 Mpc. 


5.4 The survey geometry 


The estimator presented in only provides an unbiased 
prediction of the true underlying power spectrum with¬ 
out any survey geometry effects. At intermediate and large 
scales, the measurement is affected by the shape of the 
survey, especially for high order multipoles. Given a theo¬ 
retical anisotropic power spectrum P**'®°'(k'), the observed 
power spectrum due to the effects of the survey is windowed 
through the following expression. 


■« = / 


(27r)3 

where W is defined as. 




= 7172 / d^rUsiOe' 


( 10 ) 


( 11 ) 


We refer to \W\‘^ as the window function, which satisfy 
the normalization condition, J d^k'/(27r)^ |iy(k')p = 1, im¬ 
posed by the definition of the factor I 2 in Eq. The func- 


10 


is a convolution. 


tional provided by 
Therefore, the convolution theorem can be apply making 
use of EFT techniques, which allows the computation in a 
minutes-time scale per model 

We assume that the monopole and quadrupole provide 
all the information about the full /r-shape of the power spec¬ 
trum. Thus we can write, 


ptheo.(k) ^ p (12) 


We can define a windowed power spectrum Gmultipole as, 
p(^)"'‘“ (fc) = I dMkTr(/rk)P'*‘" (k), (13) 


where fiii = ■ k. From Eq. 13 it is clear to see that both 


p(o)"'‘"- g^jjd P^^^ *1)1, have contributions from both P^°^ 


and pG)*''^-. Eq. 


13 


provides a full description of the effect 


of the window in the monopole and quadrupole. However, 
for practical reasons, is convenient to use this equation to 
calibrate a matrix that is able to relate the pG)"‘“' a giv en 


fc-bin from an arbitrary shape of P^^^ 
write the matrix elements, 

2e' -h 1 
2 

3{r)theo. 


’. Using Eq. 


13 


W“ = 


d^jL^COuk) J d3k'|W(k-k')|' 


{k')Ce{ii')eTii{ki - fc')] /P‘^'^“'“'(fci), 


(14) 


where 0th (fci — k') is a top hat function around the fci-bin: 
0TH(fci — k') is 1 when k' belongs to fci-bin and 0 otherwise. 
In order to form a window-matrix which is able to mimic 
the behaviour described by Eq. we have used 1000 k- 
bins between kf and k = 0.5 /iMpc“^ as an input of kj and 
60 output fci-bins which coincides with the fc-values where 
we measure the dataset. Using the calibrated values of the 
window matrix described by Eq. we write the windowed 
monopole and quadrupole power spectra as. 



j 

j 






(15) 
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Note that Eq. and describe the same survey window 
effect, but the latter one is much faster to be applied to 
minimization and mcmc algorithms. In this paper we al¬ 
ways use combined window that we obtain by weighting the 
individual windows of NGC and SCG by their area. Since 
this is what we do with the power spectrum measurements, 
the combined window reproduce by definition the combined 
power spectrum. 


6 PARAMETER ESTIMATION 

In this section we describe how the parameters of interest 
and their errors, including the AP parameters and /ug are 
estimated. We also present a systematic test on the power 
spectrum model presented in fusing the galaxy mocks. 

6.1 Covariance Matrices 


The matrix terms in Eq. depend, in principle, on 
the choice of theoretical power spectrum, . However 

this dependence is expected to be weak, as the term P^^^ 
appears in both denominator and numerator, and the part 
of the signal in Wij coming from P^^^ is cancelled. In 
order to quantify this, we have built two window matrices 
with different power spectra, whose values are within 2a of 
the data and no significant differences were observed. 


Fig. shows the effect of the LOWZ (left panel) and 
GMASS (right panel) DR12 samples on the power spec¬ 
trum monopole and quadrupole. The upper panel shows 
a toy model input power spectrum monopole (blue solid 
lines) and quadrupole (red solid lines) and its output ac¬ 
cording to Eq. 15 (blue and red dashed lines). The lower 
panels show the ratio between the input and the output 
power spectrum multipole. We see that for GMASS the ef¬ 
fect, both in the monopole and in the quadrupole, is smaller 
than for the LOWZ. This is due to the size of the sam¬ 
ple. The larger the sample is, the smaller is the effect of 
the mask. At fc ~ 0.2/iMpc“^, the effect of the window 
is near 0% for the monopole in GMASS, ~ 5% for the 
quadrupole in GMASS, ~ 7% for the monopole in LOWZ 
and ~ 10% for the quadrupole in LOWZ. At large scales, 
the windowed power spectrum predicted for the quadrupole 
has oscillatory behaviour, a consequence of the lack of /r- 
modes in the fc-bins at these scales; since the fe-binning 
is logarithmic, the number of fundamental fc-modes con¬ 
tained at large scales is much smaller than that at small 
scales. As a consequence, at these large scales the con¬ 
straint J dij.£.e{fi)jCei (jj.) = is not satisfied for 1,1' > 0 
with sufficient precision. The behaviour can be corrected 
by renormalising, both data and window by the sum of 
squared Legendre polynomials over the p-modes: /r-modes 
X];j-modes( m)- However, since for this paper we only use 
quadrupole data for k > 0.04/iMpc“'^, where the /r-modes 
contained in the fc-bins are large enough, the correction is 
not significant enough to be worth including. In any case, 
since the correction would have to be applied to both data 
and window, the relative ratio between theory and data 
would be unaffected by this correction. 


In this paper we do not correct the power spectrum 
multipoles for the integral constrain, which only produce a 
significant effect at scales comparable to the size of the sur¬ 
vey (see e.g. Peacock fc Nicholson|199~ I, which in our case is 
fcf ~ 0.002 /iMpc“^. Since the largest scales we consider are 
10 and 20 times smaller for the monopole and quadrupole, 
respectively, in our analysis the integral constrain is a sub¬ 
dominant component compared to other effects such as sys¬ 
tematic weights. 


The covariance matrix of the monopole and quadrupole is 
computed using the different realizations of the two sets of 
gtalaxy mocks described in §2.2[ We take into account the 
covariance of the monopole and quadrupole fc-bins and also 
the cross-covariance between these two. Each element of the 
covariance matrix is calculated from the mocks. 


m=l 

X [pf>(fc,)-(p(^'>(fc,))], (16) 

where {P^^\k)) = Prn {k)/Nm is the mean of the l- 

multipoles among realizations, and Nm is the number of 
independent realizations. The full covariance matrix for the 
monopole and quadrupole can be written in terms of the 
matrices , for £ = 0, 2, as. 


fc°° 

(^^20 


C22) ■ 


(17) 


For the QPM mocks the number of independent realization 
is Nm = 1000, whereas for MD-Patchy is Nm = 2048. In 
both cases the number of elements is much larger than the 
total number of bins, ni,, which for this work is Uf, = 120. 

Since the covariance matrix C is estimated from a set of 
mocks, its inverse is biased due to the limited number 
of realizations. We account for this effect by applying the 
correction proposed by Hartlap, Simon & Schneider (20071. 
In addition to this scaling, we have to propagate the error in 
the covariance matrix to the error on the estimated param¬ 
eters. We do this by scaling the variance for each parameter 
by the factor of eq. 18 of Percival et al. (20141. However, we 
observe that the correction due to this effect is subdominant, 
namely ^2%. 

The middle and top panels of Fig. display the cor¬ 
relation coefficient matrices, rij = Ci,j/[Ci,iCj,jY^^ for 
LOWZ and GMASS samples using the Qpm mocks and the 
MD-Patchy mocks as labeled. In all cases NGC and SGC 
have been combined into a single sample, as described in 
0 We observe that the off-diagonal terms of the auto¬ 
covariance (the covariance between monopole-monopole and 
quadrupole-quadrupole) are significantly correlated at large 
scales because of the effect of the survey geometry, for both 
monopole and quadrupole. The off-diagonal terms of the 
cross-covariance between monopole and quadrupole present 
significantly smaller correlation. 

As we go to smaller scales the auto- and cross-covariance 
off-diagonal terms are reduced for both LOWZ and GMASS 
samples because the effect of the survey window is less im¬ 
portant. For fc ^ 0.2 the auto-covariance off-diagonal el¬ 
ements start growing again because of the effect of mode 
coupling, which becomes more important as we go to smaller 
scales. At the same scales, the off-diagonal terms of cross¬ 
covariance stay very close to 0. This suggests that the mode 
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Figure 2. Top sub-panels: Effect of the window function on the monopole (blue lines) and on the quadrupole (red lines) for the LOWZ- 
DR12 sample (le ft p anel) and for the CMASS-DR12 sample (right panel). The solid lines correspond to a toy-model for ptheo.(O) 
ptheo.(2) jjj Eg. I 12 I The dashed lines correspond to and p(2)'^‘"' Eg. |l5| Lower sub-panels: Relative deviation between 

convolution with the window mask. 


coupling induces a strong correlation between close fc-modes 
in the monopole and quadrupole, but not a correlation be¬ 
tween these two statistics until k = 0.3/iMpc“'^. 

The bottom panels of Fig.|^show the ratio between the 
QPM and MD-Patchy covariances for LOWZ and CMASS 
samples, as labeled, in order to stress their differences. 
We observe that most of the off-diagonal signal is very 
noisy, although some differences can be seen for those off- 
diagonal terms close to the diagonal. For both LOWZ and 
CMASS samples, both qpm and MD-Patchy mocks pre¬ 
dict the same degree of correlation on large scales. As we 
go to smaller scales some differences arise. At scales where 
k ^ 0 . 1 , ZiMpc”^, the off-diagonal terms of the monopole 
auto-correlation matrix, tend to be more correlated in the 
QPM mocks than in the MD-Patchy mocks, while the auto¬ 
correlation matrix estimates for the quadrupole are similar. 

In Fig.j^we show the percentile diagonal error for power 
spectrum multipoles, corresponding to the fc-bins used in 
!|^ relative the power spectrum multipole amplitude of the 
data. At large scales both qpm and MD-Patchy predic¬ 
tion agree well for both the LOWZ and CMASS samples 
for both monopole and quadrupole. At small scales qpm 
mocks predict higher mode coupling in the monopole than 
MD-Patchy and therefore the errors of the power spec¬ 
trum monopole saturate before that of the MD-Patchy 
predictions. The precision of the monopole exceeds 1% for 
k ^ 0.13hMpc“^. 

For the quadrupole the relative error predicted by QPM 
and MD-Patchy mocks is very similar in all the range of 
scales studied. 

Overall, the covariances extracted from QPM and MD- 
Patchy mocks are similar and the main discrepancies are for 
the monopole at small scales. We are not able to discriminate 
which of these to sets of mocks is best. 



Figure 4. Percentile diagonal errors corresponding to the fc-bins 
used in j|^in which the monopole and quadrupole have been mea¬ 
sured for the mocks. Solid lines display the CMASS-DR12 statis¬ 
tics and dashed lines the LOWZ-DR12 ones. Red lines are the 
predictions inferred from the MD-Patchy mocks, whereas the 
blue lines are according to QPM mocks. 

6.2 Best-fit and error estimation 

We model the amplitude and the shape of the power spec¬ 
trum monopole and quadrupole through a set of 8 free pa¬ 
rameters T” = { 6 i, 62 ,/,o-g, Anoise.crpoG.aipax}, which we 
briefly describe below. 

(i) The galaxy bias is modelled using two bias parameters, 
bi and 62 as described in §5.1| The value for the non-local 
bias parameters is set by the value of bi under the assump¬ 
tion of local Lagrangian bias. 

(ii) The logarithmic growth factor /. This parameter can 
be predicted for a specific cosmological model (when the 
Qrn value is known) if we assume a theory of gravity. In this 
paper we consider this as free parameter in order to test 
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Figure 3. Correlation coefficients of the monopole-quadrupole covariance matrix from LOWZ-DR12 sample (left panels) and from the 
CMASS-DR12 sample (right panels), extracted from 1000 realizations of the QPM mocks and from 2048 realizations of the MD-Patchy 
mocks. The top panels show the results for QPM mocks, the middle panels for the MD-Patchy mocks, and the bottom panels their ratio. 


potential deviations from GR, or, if assuming GR, for not (iv) The amplitude of primordial dark matter power spec- 

using a prior on the value. trum, ng. 


(iii) The AP parameters, ay and ax- As described in 1 5.3 


through varying ay and a± we are able to parametrise the 
anisotropy generated in the power spectrum multipoles by 
assuming an incorrect value of flm- 


(v) The amplitude of shot noise, Anoise, as described in 


(vi) The Fingers-of-God parameter, ctfoG, introduced in 
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Note that although we allow / and as to vary indepen¬ 
dently, these two parameters are highly degenerate when 
constrained only using the power spectrum multipole^ We 
will only quote the measured combined quantity /erg. Simi¬ 
larly, we will only report bias, and b 2 as. 

The other cosmological parameters, such as Qrn, the 
spectral index and the Hubble parameter h are fixed at 
the fiducial values described in §2.3| during the fitting pro¬ 
cess. In we will vary the value of fim to analyse the 
effect of such a change can produce in fas ■ In this paper we 
always perform the parameter fitting process to the com¬ 
bined NGC-I-SGC sample, both for the mocks and for the 
data. 

In order to preform the parameter estimation, we as¬ 
sume the monopole and quadrupole are drawn from a mul¬ 
tivariate Gaussian distribution and use, 

X^(T-) = [AP(T-)][C-i][AP(T-)l‘ (18) 

where, AP(T’) is the vector whose elements contain the dif¬ 
ference between the data and the model for the power spec¬ 
trum monopole and quadrupole, and C~^ is the inverse co- 
variance matrix. By minimizing the function respect to 
T- we obtain the best-fit set of parameters. The errors as¬ 
sociated to each parameter are computed by exploring the 
likelihood surface when a specific parameter is fixed and the 
other parameters can vary freely. The likelihood surface is 
explored using a simplex minimization algorithm. In order 
to ensure that the minima found are global and not local, 
we run the algorithm multiple time with different starting 
points and different starting variation ranges, and check for 
convergence. 

6.3 Tests for the galaxy mocks 

In this section we test for potential systematics of the model 
presented in using the MD-Patchy galaxy mocks. We 
start by comparing the power spectrum multipoles for the 
mean of the 1000 and 2048 realizations of the QPM and MD- 
Patchy mocks, respectively, with the measurements of the 
DR12 data. This is shown in Fig. where the power spec¬ 
trum monopole and quadrupole are divided by a smooth 
linear power spectrum. The left and right panels show the 
measurements corresponding to LOWZ and CMASS sam¬ 
ples, respectively. Blue symbols represent the measurements 
for the data corresponding to the monopole, whereas the red 
symbols represent the quadrupole. The black dashed and 
dotted lines show the mean values of the mocks for QPM and 
MD-Patchy, respectively. 

For the LOWZ sample the MD-Patchy mocks de¬ 
scribe accurately the data measurements for both monopole 
and quadrupole. The QPM mocks show some discrepancies 
with the data: QPM mocks systematically overestimate the 
monopole data by ~ 10% at k ^ 0.02, /iMpc“^ and fail to de¬ 
scribe the quadrupole for fc ^ 0.1 /iMpc“^. For the CMASS 

® / and (Tg are fully degenerated only in the large-scale limit, 
when Kaiser formula is still valid. At smaller scales the non-linear 
corrections break this degeneration. However, the signal-to-noise 
of the power spectrum at these scales is not sufficiently high to 
make reasonable constrains on / and erg alone, since the degen¬ 
eration is poorly broken. Because of this we refer to these two 
parameters as highly degenerate but not fully degenerate. 


sample, MD-Patchy mocks describe well the monopole 
and quadruple data. QPM mocks tend to overestimate both 
statistics by few percent. The version of QPM mocks used in 
this paper was not designed to describe redshift space dis¬ 
tortion features with few percent accuracy, as this version of 
the MD-Patchy mocks was. Because of this, we will only 
use the MD-Patchy mocks multipole measurement to test 
the modelling of RSD. 

We combine the model described in with the mea¬ 
sured mock power spectrum monopole and quadrupole, av¬ 
eraged over the 2048 realization of the MD-Patchy mocks, 
in order to recover the /erg parameter. Since we know the 
input cosmological parameters for the mock simulations, we 
can compare the obtained value with the expected one, and 
thus test which is the precision of the model when recov¬ 
ering fas- Note that the measurements of the monopole 
and quadrupole of the mocks have been performed using 
the fiducial cosmology, which is different to the cosmology 
of the mocks as described in |2.2[ 

Figj^displays in blue symbols linked by solid blue lines 
the obtained fas, and a± parameters as a function of the 
minimum scaled used to fit the model to the measurements 
of the MD-Patchy mocks. The error-bars correspond to the 
data error-bars scaled by a volume factor of V2048 in order 
to account for the volume difference. In order to mimic the 
data analysis, the large scales cuts fcmin = 0. 02 hMpc ^ and 
0.04 /iMpc“^ for the monopole and quadrupole, respectively, 
have been applied. In addition to Q||, ax, / and erg, we al¬ 
low hi, & 2 , Anoise and (JfoG to vary, as we do for the data. 
The expected values for fas are shown in black dashed lines. 
For reference, 1% and 3% deviations from the expected fas 
are also shown in dot-dashed and dotted black lines, respec¬ 
tively. 

For the LOWZ sample we see that the recovered value 
of fas agrees to ^ 2% accuracy with the true one for 
fcmax ^ 0.16 hMpc“^. We observe that if the truncation scale 
is fcmax = 0.15 hMpc“^, the model produces a systematically 
low value of fas. However, this does not occur when the 
truncation scale is higher than this value. In particular, if 
the truncation scale is fcmax 0.20/iMpc“'^ the accuracy of 
the obtained fas is < 1%, which is around 10 times smaller 
than the expected statistical error for this sample. We ob¬ 
serve that ax is constrained with < 1% accuracy, at all 
studied scales, ay presents ~ 3% deviation from its fiducial 
value at fcmax = 0.24hMpc“^, and decreases to ^ 1% as the 
scale of truncation is increased. 

For the CMASS sample we see that the recovered value 
of fas agree with < 1% accuracy for truncation scales within 
the whole range of scales studied. Given the statistical errors 
for the CMASS sample, this systematic error is around 8 
times smaller. As for the LOWZ case, ax is recovered within 
Si 1% accuracy and ay within ^ 2%, in all the range of 
studied scales. 

Since the observed systematic errors for both LOWZ 
and CMASS samples are much smaller than the statistical 
errors obtained for the data, we do not correct fas by any 
systematic shift. 

Fig. in shows the distribution of parameters f{z)as{z), 
H[z)rs{zd) and DA{z)/rs{zd) for the 2048 realizations of 
the MD-Patchy mocks in the LOWZ and CMASS sam¬ 
ples. The blue crosses represent the best-fit parameter set 
for each realization of the mocks, whereas the red cross is 
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Figure 5. Performance of the QPM (black dashed lines) and MD-Patchy (black dotted lines) mocks compared to the data, blue squares 
for the monopole and red circles for the quadrupole. The left and right panels display the results for the LOWZ and CMASS samples, 
respectively. For clarity, the amplitude of the power spectrum multipoles have been normalised by a smoothed linear power spectrum, 

-^smooth • 



0.15 0.16 0.17 0.18 0.19 0.20 0.21 0.22 0.23 0.24 
kmaxIhMpc'l 



0.15 0.16 0.17 0.18 0.19 0.20 0.21 0.22 0.23 0.24 

kmax[HMpC-'] 


Figure 6. Best-fit qx, 0|| and fag parameters as a function of the minimum scale (maximum k) used for describing the LOWZ- 
(left panels) and CMASS (right panels) monopole and quadrupole power spectra predicted by the average of 2048 realizations of the 
MD-Patchy mocks (blue squares linked by solid lines). In order to mimic the analysis of the data, the largest scales used for the fit are 
^min = 0.02/iMpc“^ for the monopole and fcmin = 0.04/iMpc“^ for the quadrupole. All values normalised by the corresponding true 
expected value. Horizontal dashed, dot-dashed and dotted lines show the 0%, 1% and 3% deviations, respectively, with respect to the 
corresponding fiducial value. 


for the data. The shapes of the clouds formed by the blue 
crosses provide information of how these parameters are de¬ 
generate. Thus, these plots show that the measurements of 
fag, H{z)rs{zd) and DA{z)/rs{zd) are strongly correlated. 
Assuming that the data results can be approximated by a 
multivariate Gaussian likelihood, in §7.5| we provide the co- 
variance matrix among these 3 parameters for the measure¬ 
ments from the data. 

Fig-Hl shows the two-dimensional parameter space of 
b\(jg and fag for the 2048 realizations of the MD-Patchy 
mocks in the LOWZ and CMASS samples. 

We see that when the AP parameters are set to 1, the 
scatter between best-fit values from different mocks reduces 
significantly forming a small cloud (green crosses) which is 


contained by the blue one. The red crosses represent the 
actual values for the data when the MD-Patchy covariance 
is used. In both cases we see that the data points lies well 
within the clouds formed by the different realizations of the 
MD-Patchy mocks. 


7 BOSS DR12 MEASUREMENTS 

The values of the parameters of the model corresponding 
to the fits to the data presented in Fig. are listed in Ta¬ 
ble where the minimum cut-off scale used for the fit is 
fcmax = 0.24hMpc“^. We observe that the difference be¬ 
tween parameters (using QPM and MD-Patchy covariances) 
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Figure 7. Two dimensional distribution of f{z)as{z), H{z)rs{zd) and Da{z.)/ rs{zd)■, for the LOWZ and CMASS samples, left and right 
panels, respectively. The blue points represent the best-fit solutions of 2048 independent MD-Patchy mocks realizations. The red crosses 
are the results for the data. The scatter of the mocks gives an idea of how these parameters are degenerate. Both mocks and data have 
been analysed using the MD-Patchy covariance. 


□ R12 LOWZ MD-Patchy DR12 CMASS MD-Patchy 





Mocks (AP test) 




X Mocks (no AP test) 

0.8 


0.8 
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^ 0.6 
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^ 0.4 
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Figure 8. Two-dimensional distribution for the parameter space corresponding to bias and fag. Left and right panels display the 
results for the LOWZ and CMASS samples. Blue x-symbols correspond to the best-fit solution for each realization (out of 2048) of 
the MD-Patchy mocks when all parameters are freely varied. The green -|~symbols correspond to the best-fit solution when the AP 
parameters, 0 || and ay are set their fiducial values, respectively. The red X-symbol display the best-fit solution for the data sample 
when all the parameters are varied. The red -|—symbol display the best-fit solution for the data sample when the AP parameters have 
been set to their fiducial values. In all the cases the covariance matrix have been inferred using the MD-Patchy mocks. 


and their error-bars is not significant. This indicates that the 
impact of the different covariance matrices in the parame¬ 
ter estimation is sub-dominant compared to other effects, 
such as the scale where the fitting process is truncated, or 
systematics of the model itself. 

Table displays the best-fit results at fcmax = 
0.24/iMpc“ , with H{z)rs{zd) and DA{z)/rs{zd) set to the 
fiducial cosmology prediction. As in Table the differences 
between the parameters and their errors obtained using ei¬ 
ther QPM and MD-Patchy mocks are not significant. 

In Fig. 1^ we show how the best-fit values of /erg and 
the AP parameters, ay and ay, change with the minimum 
scale considered. 

From the left panels of Fig. we observe that when 


varying fcmax from 0.15/iMpc“^ to 0.24hMpc“^, the best- 
fit value of /(jg is very stable compared to the size of the 
error-bars. This behaviour applies when the AP parameters 
are varied and also when they are fixed to their fiducial 
value When the AP parameters are varied the best-fit value 
of /(Tg is consistently low by ~ Icr respect to Planckl5+GYi 
prediction. We also observe that, whereas ay is consistent 
with unity for all values of fcmax studied, ay presents a low 
value of 2 ct — 3cr respect to the fiducial cosmology predic¬ 
tion. When the AP parameters, ay and Qy, are set 1, the 
LOWZ-sample /erg parameter, is shifted ~ Icr upwards and 
is more in accordance with the PlancklS cosmology predic¬ 
tion. We observe that this effect is independent of the model 
truncation scale. In §7.1|we will reanalyze the data assuming 
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Figure 9. Power spectrum monopole and quadrupole best-fit parameters as a function of the minimum scale considered, fcmax- The 
results for the LOWZ- (left panels) and CMASS-DR12 (right panel) BOSS data are shown in red solid lines in combination with red 
filled symbols. In the fas panel we also show in dashed blue lines and empty blue symbols the best-fit values when the AP parameters 
have been set to their fiducial value. In this case, the symbols have been displaced horizontally for visualization reasons. For all the cases 
the largest scales used for the fitting are kmin = 0.02/iMpc“^ for the monopole and fcmin = 0.04hMpc“^ for the quadrupole. In the fas 
sub-panels we show as a black solid line the predictions from Planckld. The cyan and magenta bands represent the la and 2a error-bars, 
respectively, around the Planckld best-fit solution. 


Sample (Cov.) 

LOWZ (qpm) 

LOWZ (MD-Patohy) 

CMASS (qpm) 

CMASS (MD-Patchy) 


0.392 ±0.061 

0.395 ± 0.064 

0.445 ± 0.038 

0.442 ± 0.037 

H{z^fi)rs{zd) [lO^kms”^] 

11.48 ±0.55 

11.33 ±0.56 

13.99 ±0.44 

13.84 ± 0.43 

DA{z)/rs{zd) 

6.38 ±0.18 

6.33 ±0.19 

9.43 ±0.15 

9.42 ± 0.15 


1.025 ±0.052 

1.039 ±0.054 

0.977 ± 0.030 

0.987 ± 0.030 


0.945 ± 0.027 

0.938 ± 0.028 

0.999 ±0.016 

0.998 ±0.016 

bie^siZeff) 

1.283 ±0.032 

1.279 ±0.037 

1.218 ±0.022 

1.225 ± 0.020 

b2(^s{Zef!) 

-0.19 ±0.64 

-0.38 ±0.36 

0.67 ± 0.74 

0.40 ± 0.66 

A 

■^noise 

-0.30 ±0.22 

-0.43 ±0.21 

-0.041 ± 0.078 

-0.057 ± 0.093 

0'FoG(2;eff)[Mpc/l-l] 

3.94 ± 0.56 

4.23 ±0.56 

3.35 ±0.32 

3.42 ± 0.31 

x^/d.o.f 

29.62/(53-8) 

31.48/(53-8) 

26.168/(48-8) 

33.661/(48-8) 


Table 1. Best-fit parameters obtained from fitting the monopole and quadrupole BOSS DR12 data using the theoretical approach 
described in ^ The two first columns are the parameters obtained from fitting the LOWZ-DR12 data, whereas the third and forth 
columns are obtained from fitting the CMASS-DR12 data. The first and third columns are the parameters obtained when the covariance 
matrix is inferred from QPM mocks, whereas the second and forth columns parameters are obtained when the covariance matrix is 
inferred from MD-Patchy mocks. For all cases, the minimum scale used is fcmax = 0.24/iMpc“^ and the largest scales used are 
^min = 0.02hMpc“^ for the monopole and fcmin = 0.04hMpc”^ for the quadrupole. The error-bars represent la deviations and have 
been inferred from the analysis of the likelihood function as it is described in ^6.2| The cosmological parameters, fas, II(z)rs( 2 :^) and 
Dj[(z)/rs(z^) are correlated, so we encourage using the multivariate Gaussian likelihood presented in ^7.5| 


Sample (Cov.) 

LOWZ (qpm) 

LOWZ (MD-Patchy) 

CMASS (qpm) 

CMASS (MD-Patchy) 

/(^eff)o'8(^efr) 

fel(T8(2eff) 

b2crs(zes) 

A 

■^noise 

o-FoG[Mpch“l] 

0.476 ± 0.043 
1.290 ± 0.016 
-0.23 ±0.95 
-0.32 ±0.22 
4.04 ± 0.61 

0.494 ± 0.045 

1.283 ±0.017 
-0.37 ± 0.60 
-0.41 ± 0.21 

4.27 ± 0.59 

0.434 ± 0.023 
1.236 ±0.012 

0.38 ±0.60 
-0.044 ± 0.089 
3.41 ± 0.31 

0.438 ±0.021 

1.236 ±0.012 

0.25 ±0.56 
-0.064 ± 0.086 

3.45 ±0.30 

X^/d.o.f 

33.45/(53-6) 

36.08/(53-6) 

26.96/(48-6) 

33.99/(48-6) 


Table 2. Best-fit parameters obtained from fitting the monopole and quadrupole BOSS DR12 data as it has been described in Tablef^ 
where in this case the AP parameters, D/^{z)Its{^z^ and Hi^z)rs{Z(i) have been set to their fiducial value, 77(-2iowz)'^(^d) = 11.773 kms“~ 
DA{z\o-^^)lrs{zd) = 6.7466, H{zcuiass)r{zd) = 13.663 ■ 10®kms“\ Da(z\o^^,)/ ra{zd) = 9.4418. 
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different fiducial models for and we will discuss how this 
tension changes. So far we can say that the LOWZ sample 
data has a preference for low values of fag and q_l (which is 
translated into low values of DA{z)/rs{zd), according to the 
definitions in (5.31, which is within 1 and 2a of the fiducial 


values when GR is assumed. Both parameters are correlated 
as we will show in §7.5| so when qx is fixed to its fiducial 
value, fag is automatically shifted towards a solution which, 
in this case, is consistent with Planckl5+GK within \a. 

From the right panels of Fig. we see that the values of 
AP parameters are both consistent with their fiducial val¬ 
ues within la for all fcmax, suggesting that there is no strong 
tension between the fiducial cosmology and the actual. Also, 
the inferred fag value when the full AP test is performed, 
is consistent with PlancklS+GK within la. By setting the 
AP parameters to their fiducial values, we do not observe 
any significant change on the fag parameter, besides a re¬ 
duction of the error-bars, which increase the tension with 
PlancklS+GP by up to 2a. 

In Tablej^we present the final measurements of this pa¬ 
per. They have been obtained by averaging the results pre¬ 
sented in (Q(Table[^and|^ using the QPM and MD-Patchy 
covariances. We show the results for both the LOWZ and 
CMASS sample where the NGC and SGC have been com¬ 
bined into a single measurement. The main result of this pa¬ 
per is /(0.57)(T8(0.57) = 0.448 ± 0.038 for the GMASS and 
/(0.32)(js(0.32) = 0.402T0.060 for the LOWZ sample when 
the AP test is performed; and /(0.57)crg(0.57) = 0.438 ± 
0.022 for the CMASS and /(0.32)cr8(0.32) = 0.497 ± 0.0436 
for the LOWZ sample then the AP parameters have been 
tuned to their fiducial value . In the following subsections 
we analyse how these results are affected by a change in the 


cosmological models (( 7.11, how these results compares with 


other values of f{z)ag{z) extracted from papers based on the 
Data Release 11 of the BOSS sample ((7.31, and how these 
results compare to other values of f{z)ag{z) obtained from 
other redshift surveys ( |7.4[ |. We also present the multivari¬ 
ate likelihood surface at fag, H{z)rs{zd) and DA{z)/rs{zd) 


((7.51. 


T.l Dependence with cosmology 

In this section we study how sensitive the measurements 
of a^, Oil and fag are respect to the assumed cosmology 
for converting redshifts into comoving distances. In order to 
test this effect we have measured the galaxy power spec¬ 
trum multipoles of the data assuming two additional cos¬ 
mological models, one with Qm ~ 0.292 and another with 
= 0.332 (see table 4 of Gil-Marfn et al.|2015a for details 
about these two extra cosmologies, listed as H-Planckl3 and 
L-Planckl3). These two values of tlm are at Icr tension of the 


best-fit value reported by Planck Collaboration et al. (20141. 
Our power spectrum models are as described in but based 
on the linear power spectrum of the assumed cosmological 
model (see fig. 6 in Gil-Marfn et al.|2015a for a comparison 
between these linear dark matter power spectra). In order to 
be consistent, we have re-analysed the MD-Patchy mocks 
assuming the new value of Dm for converting redshifts into 
comoving distances, in order to recompute a new covari¬ 
ance matrix, consistent with the new measurement of the 
monopole and quadrupole. In Fig. we show the values 
of a_L, a II and fag obtained as a function of the chosen 


value for Dm for the LOWZ and CMASS sample in the left 
and right panels, respectively, using the corresponding MD- 
Patchy covariance. The red, blue and green points show the 
results of assuming Dm = 0.292, 0.310, 0.332, respectively. 
The filled symbols represent the results when the full AP 
test is performed, whereas the empty symbols represent the 
values fag obtained when the AP parameters have been set 
to their fiducial value (in this case Q|| and a±^ are set to 
1). In the bottom sub-panel panel, the dotted horizontal 
coloured lines show the corresponding value of fag when 
GR is assumed and the value of <78 is set to 0.815 (which 
is the fiducial value for the PlancklS cosmology). In the a±_ 
and ay panels, the horizontal black dotted line represent the 
fiducial value = ay = 1. 

For the CMASS sample, the results do not strongly de¬ 
pend on the fiducial cosmology chosen. For both AP and 
non-AP cases, the changes in fag as a function of the cos¬ 
mological model assumed are typically of ~ 0.5a. Further¬ 
more, when the Dm = 0.31 model is assumed, the tension 
between the observed fag and its Planckl5+GK prediction 
minimises. Therefore, CMASS data suggests that when GR 
is assumed as a theory of gravity. Dm = 0.31 is the most 
likely model among those studied. Also, for this case the 
AP parameters are the closest to the fiducial values. 

For the LOWZ sample, when the AP test is performed 
we observe a ~ 2 (t tension for D))f = 0.293 and 0.310 on 
fag. The tension is even higher when D^ = 0.332. Sim¬ 
ilar findings apply to the best-fit values of a±_ parameter, 
whereas the best-fit value for ay is similar for the 3 studied 
models. When the AP test is turned off, and consequently 
ay and ax are set to their fiducial values, we observe that 
the fag tension is reduced to well within Icr for all models 
studied. This sug gests, that the cause of the fag tension ob¬ 
served in Fig. [^cannot be explained by the choice of D))f. 
Thus, according to LOWZ data, the observed DA{z)/rs{zd) 
parameter is in ^ 2 ct tension with its fiducial value, which 
induces a ^ 2cr tension in the observed fag compared to its 
Planckl5-\-GYl prediction, due to the correlation between pa¬ 
rameters. This tension is relaxed when DA{z)/rs{zd) is set 
to its fiducial value. In any case, the Dm = 0.31 model is the 
choice that presents less tension between the observed fag 
and its Planckl5+GK prediction. 

Overall, we conclude that both CMASS and LOWZ 
galaxies data is in agreement with the fiducial model D))f = 
0.31 consistent with PlancklS data: for the CMASS sample 
the tension is below la and for the LOWZ sample within 2cr. 
For both samples, D^jf = 0.31 is the model that minimises 
the tension between the observed fag and the predicted by 
PlancklS+GP. 


7.2 Comparison of fag with other DRll BOSS 
Data Releases analyses 

In this section we compare our results on fag with other 
studies of RSD based on the previous Data Release 11 
(DRll) of the BOSS sample. DRll contains ~ 10% less 
galaxies than DR12, however, we expect that this change is 
not sufficiently large for producing signifiant changes (other 
than systematics) in the obtained best-fit of fag. In Fig. 
f^we show the measurements of fag based on the follow¬ 


ing DRll works: (Chuang et al. 2013 Beutler et al. 2014 
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LOWZ 

CMASS 

LOWZ (no-AP) 

CMASS (no-AP) 

/( 2 eff)'rg(Zeff) 

0.394 ± 0.062 

0.444 ± 0.038 

0.485 ± 0.044 

0.436 ± 0.022 

J 7 ( 2 eff)rs(zd) [ 10 ®kms“l] 

11.41 ± 0.56 

13.92 ±0.44 

11.773 

13.663 

DA{z)/rs{za) 

6.35 ±0.19 

9.42 ±0.15 

6.7466 

9.4418 

“||( 2 eff) 

1.032 ±0.053 

0.982 ±0.031 

1 

1 

Oj_( 2 eff) 

0.942 ± 0.027 

0.998 ±0.016 

1 

1 

biagiz^s) 

1.281 ± 0.035 

1.222 ± 0.021 

1.287 ± 0.017 

1.236 ±0.012 

b2ag(zgff) 

-0.29 ±0.50 

0.53 ±0.70 

-0.30 ±0.78 

0.32 ± 0.58 

4 

■^noise 

-0.36 ±0.22 

-0.049 ±0.086 

-0.36 ±0.22 

-0.054 ± 0.077 

o-FoG[Mpc/l“l] 

4.08 ±0.56 

3.39 ± 0.32 

4.16 ±0.60 

3.43 ±0.31 


Table 3. Combined best-fit parameters for the LOWZ- and CMASS-DR12 samples. The first column display the LOWZ results; the 
second column display the CMASS results; in both cases the AP test is performed. The third and forth column display the results for 
LOWZ and CMASS, respectively, when no-AP test is performed. The combination has been performed taking the average of the best-fit 
values of parameters obtained when the covariance matrix is extracted from either QPM and MD-Patchy mocks, as it has been listed in 
Tables (when AP test is performed) and(when no-AP test is performed). The displayed errors are also taken as the average as the 
ones listed in TablesandThe minimum scale for the fit is fcmax = 0.24hMpc“^, and the large scale cuts are fcmin = 0.02/iMpc“^ 
for the monopole and fcmin = 0.04hMpc“^ for the quadrupole. 



0.293 0.31 0.332 



0.293 0.31 0.332 


n 


m 


n 
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Figure 10. The left and right panels show the best-fit ax (top sub-panel), ay (middle sub-panel) and fag (bottom sub-panel) value for 
LOWZ and CMASS sample, respectively, as a function of the chosen cosmological model to convert redshift into distances. The different 
colours show 3 different models, Qrn = 0.293 (red), Qrn = 0.31 (blue) and Qrn = 0.332 (green). The filled squares show the results when 
the AP parameters are varied, whereas for the empty squares the AP parameters have been set to their fiducial values. In the bottom 
sub-panel, the horizontal dashed lines show the fag prediction assuming GR (/(z) = , with 7 = 0.545) and ag = 0.815. 


Samushia et al.||20141 |Sanchez et aL|2014 

|Reid et aL||2014[ 

Alam et al.||2015b Gil-Marm et al.|2015a 

I. 


(i) Chuang et al. (20131 analyse the DRll LOWZ 


and CMASS two-point correlation function monopole 
and quadrupole in configuration space in the range 
of 56Mpc/i“^ ^ r ^ 200Mpch“^ for CMASS and 

32Mpch,“^ ^ r ^ 200Mpch“^ for LOWZ . Their work 
include fag measurements for both LOWZ and CMASS 
when the AP parameters are varied, /(0.57)(jg(0.57) = 
0.354 ±0.059 and /(0.32)crg(0.32) = 0.384 ± 0.095. 


(ii) Beutler et al.| (20141 compute the DRll CMASS 


power spectrum monopole and quadrupole in fc-space in 
the range of 0.01 hMpc“^ 5; fc 0.20hMpc“^. They re¬ 
port /(0.57)(jg(0.57) = 0.419T0.044 for the CMASS sample 
varying the AP parameters as well. 


(in) Samushia et al. (2014) compute the DRll CMASS 
two-point correlation function monopole and quadrupole 


in configuration space in the range of 24Mpc/i“^ ^ r ^ 
152Mpc/i“^. When they analyse the CMASS sample they 
obtain /(0.57)(jg(0.57) = 0.441 ± 0.044 when the AP test is 
performed, and /(0.57)(Tg(0.57) = 0.447 ± 0.028 when the 
AP parameters are tuned to their fiducial value. 


(iv) Sanchez et al. (2014) analyse the two-point correla¬ 
tion function and the clustering wedges, parallel and per¬ 
pendicular, for the DRll LOWZ and CMASS sample us¬ 
ing scales of 40Mpc/i“^ ^ r ^ 180Mpc/i“^ with AP 
test. They report /(0.57)(jg(0.57) = 0.417 ± 0.045 and 
/(0.32)crg(0.32) = 0.48 ± 0.10. 


(v) Reid et al.| (2014) perform a small scale analysis in 
the range 0.8Mpc/i“ ^ r ^ 32Mpc/i“^ using halo oc¬ 

cupation distribution and Planck cosmology. They report 
/(0.57)crg(0.57) = 0.452 ± 0.011 for the CMASS sample. 
Note that the errors on this measurements are considerably 
smaller, but at the same time, they rely in significant mod- 
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LOWZ sample CMASS sample 




Figure 11. The constrains on fag from LOWZ- (left panel) and CMASS-DRll samples (right panel) are displayed in black circles 
(|Chuang et al.|2013[ [Beutler et al.|2014[ [Samushia et al.|2014[ [Sanchez et al.|2014[ [Reid et al.|2014| [Alam et al.|2015b| [Gil-Mari'n et al?] 

|2015a| |. In red squares the predictions LOWZ and CMASS-DR12 samples, respectively, as they are listed in TableEmpty symbols are 
analysis with no-AP test, whereas filled symbols represents analysis where AP-test is performed. The cyan and magenta bands show the 
1 and 2(T-range, respectively, allowed by Planck TT+lowP-|-lensing in the base ACDM-I-GR model jPlanck Collaboration et al.|20f^. 


elling and cosmological assumptions, such as tuning the AP 
parameters to their fiducial value. 


(vi) Alam et al. (2015bI compute the DRll CMASS 
two-point correlation function monopole and quadrupole 
in configuration space in the range of 30Mpc/i“^ ^ r 
126Mpc/i~^. Their analysis includes AP results and report 
/(0.57) cr8(0.57) = 0.462 ± 0.041 for the CMASS sample. 

(vii) [Gil-Marm et al.| ( |2015a| compute the DRll CMASS 
power spectrum monopole and bispectrum monopole in the 
range of O.OShMpc”^ ^ fc ^ 0.17/iMpc“^. They report a 
value of /° "'^0.57 )ct 8(0.57) = 0.582 ± 0.084 using RSD and 
no AP effecl|j When this value is combined with jP'anckia _ 
0.777 ( where PlancklS is the fiducial cosmology used in 
that paper; see table 4 of [Gil-Marm et al.||201^ , a fag 
value can be obtained, /(0.57)(T8(0.57) = 0.504T0.069. Note 
that this measurement is not based on the power spectrum 
quadrupole and relies on different cosmological assumptions. 

The observed differences between the DRll results are 
expected to be due to different systematics in the differ¬ 
ent models, scales and statistics considered to describe the 
shape of the power spectrum or two-point correlation func¬ 


tion. With the exception of Chuang et al. (20131 in CMASS, 


all the results are consistent within 1 and 2 (t. All the results 
that perform an AP analysis have a well similar error-bars. 

For the CMASS sample, our fag measurement is con¬ 
sistent within Ict with the DRll reported values, with the 


exception of Chuang et al. (20131, which is 1.5cr below our 
findings and ~ 2a below PlancklS prediction. In particu¬ 
lar, our results are very close to those by [Samushia et al.j 
(20141 and Alam et al. {2015bI, which are also consistent 


with PlancklS data within la. When the AP parameters 
are tuned to their fiducial value, our result on fag is also 


^ The quantity constrained by the power spectrum monopole 
in combination with the bispectrum monopole is not / times ag, 
but / to the power 0.43 times ag. For a more detailed discussion 
about this topic see|Gil-Marm et al.||[^15a[ 


very close to that of Samushia et al. (2014) and consistent 
within la to the Reid et al. (2014) results. 

For the LOWZ sample our result lies between the best- 
fit values of Chuang et al. (2013) and Sanchez et al. (2014). 


7.3 Comparison of the geometrical AP 

parameters with other DR12 BOSS Data 
Releases analyses 


In this section we compare our results on the geometrical AP 
parameters, H{z)rs{zd) and DA{z)/ra{zd) with other stud¬ 
ies based on the analysis of the anisotropic 2-point function 
statistics of the Data Release 12 (DR12) of the BOSS sam¬ 
ple. 


In Table|^we compare the results from this paper (RSD- 
PS) with the BAO-based analysis of the anisotropic power 
spectrum (BAO-PS) by [Gil-Marm et al!] ( [2015b[ ) and the 
BAO analysis of the anisotropic correlation function (BAO- 
CF) by [Cuesta et al. (2015). The BAO analyses provide re¬ 
sults based on the same dataset used by the RSD analysis, 
which we refer to as pre-reconstruction analyses (Pre-BAO), 
and on the dataset obtained after applying the reconstruc¬ 
tion algorithm, which we refer as post-reconstruction anal¬ 
yses (Post-BAO). 

The results of the pre-reconstruction analyses, Pre- 
BAO-PS, Pre-BAO-CF and RSD-PS, are very consistent 
with each other for both LOWZ and CMASS samples, 
and for both H{z)rs{zd) and DA{z)/rs{zd). For LOWZ, 
H{z)rs{zd) is almost exactly the same for Pre-BAO-PS, 
Pre-BAO-CF and RSD-PS, and DA{z)/rs{zd) Pre-BAO-PS 
and RSD-PS are < 0.5 (t apart. For CMASS, P[{z)rs{zd), 
there is < 0.5cr between RSD-PS and Pre-BAO results, and 
the Da{z)/ ra{zd) result is also very consistent, being the 
same for RSD-PS and Pre-BAO-PS, and very similar with 
Pre-BAO-CF as well (< 0.5 (j). In general, the small dif¬ 
ferences observed between Pre-BAO-CF and Pre-BAO-PS 
could come from i) observational systematics (such as pho¬ 
tometric calibration, systematic and fibre collision weights) 
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Sample 

Method-Statistic 

H(z)rs{zd) 

DA(z)/ra{zd) 


Post-BAO-PS 

Post-BAO-CF 

11.64 ±0.62 

11.65 ± 0.84 

6.85 ±0.17 
6.67 ± 0.13 

LOWZ 

Pre-BAO-PS 

Pre-BAO-CF 

11.44 ±0.75 
11.44 ±0.74 

6.48 ±0.27 
6.40 ± 0.37 


RSD-PS 

11.41 ± 0.56 

6.35 ±0.19 


Post-BAO-PS 

Post-BAO-CF 

14.56 ±0.38 
14.75 ± 0.54 

9.42 ±0.13 
9.52 ±0.14 

CMASS 

Pre-BAO-PS 

Pre-BAO-CF 

14.10 ±0.65 
14.14 ±0.71 

9.42 ±0.22 
9.51 ± 0.19 


RSD-PS 

13.92 ± 0.44 

9.42 ±0.15 


Table 4. Geometrical AP parameters for LOWZ and CMASS 
samples: H{z)rs{z^) (in 10 ®kms“^ units), Dji{z)/rs{zii), in¬ 
ferred from the BAO power spectrum analysis (BAO-PS) l |Gil-| 
|Marm et al.||2015b i, BAO correlation function analysis (BAO- 
CF) ( |Cuesta et al.| 2015^ , and the RSD power spectrum analy¬ 
sis (RSD-PS) (this work). For the BAO analyses we present the 
pre- and post-reconstruction results, Pre-BAO and Post-BAO, 
respectively. Note that the RSD results are essentially a pre¬ 
reconstruction analysis, since the RSD signal is removed in the 
reconstruction process. 


that may affect differently the measurement of the corre¬ 
lation function and the power spectrum, ii) from the fact 
that the correlation function and power spectrum (in finite 
scale-ranges) do not contain exactly the same information, 
or iii) simply that the statistical noise affects the measure¬ 
ments differently in the same way as a different binning of 
the data. The differences between Pre-BAO-PS and RSD- 
PS could have their origin in systematics of the model it¬ 
self, or observational systematics that enter differently in 
the BAO and RSD modelling. In any case, these systemat¬ 
ics are smaller than O.Scr. 

On the other hand, the results between post¬ 
reconstruction analyses and pre-reconstruction (including 
RSD-PS) are more significant, but always less than 2a. Re¬ 
construction is known to reduce non-linear effects on the 
BAO signal, and some of the discrepancy could be due to 
this. This is discussed in more detail in both ICuesta et ^ 
(20151; Gil-Marin et al. (2015b I. For this reason, we recom¬ 
mend using the post-reconstruction results of DA{z)/rs{zd) 
and H(z)ra{zd), either from Post-BAO- CF, Post-BAO-PS, 
or the consensus value between both results presented in 
both papers, for constraining cosmological models. In case 
also the value of fa8 is wanted to be used (in combination 
with the geometrical AP parameters) for constraining cos¬ 
mological models, this paper is, at present, date the only 
BOSS DR12 paper that provides these measurements for 
the CMASS and LOWZ samples. 


T.4 Comparison with other galaxy surveys 

In this section we compare our measurement on fag for 
the LOWZ and CMASS samples with the reported fag val¬ 
ues of other galaxy surveys at different redshifts, and with 
PlancklS predictions. 

Fig. displays the different RSD analyses of redshift 
galaxy surveys. The analyses with filled symbols solve simul¬ 
taneously the RSD and AP effect, whereas those represented 
with empty symbols only consider RSD, keeping the cosmol¬ 
ogy fixed. 



Figure 12. Gonstraints on fag from several redshift surveys in 
the base of ACDM, with f{z) = orange circle ( 6 dFGRS 

by |Beutler et al!]||2012[|); gr ay triangle (SDSS Main Galaxy Sam- 
ple by |Howlett et al.| 2015[l); gr een inverse triangle (SDSS Lumi¬ 
nous Red Galaxies by Oka et aT)] ( |2014| ); cyan diamonds (WiggleZ 
by [Blake et al.| ||2012 i); and purple pentagon (VIPERS by jde la] 
[Torre et al.j | |2013| l). In red squares the results from BOSS-DR12 
for the LOWZ and CMASS sample according to Table Filled 
symbols represent the fag measurement when the RSD and the 
AP effect is considered, and filled symbols when only RSD effect 
is considered. For the empty red squares the z-position has been 
slightly displaced for clarity. The red, blue and green bands show 
the l(T-range allowed by Planck TT-|-lowP-|-lensing in the base 
ACDM model | [Planck Collaboration et al.[201^ when 7 = 0.420, 
7 = 0.545 (GR) and 7 = 0.680, respectively. 


(i) The GdFGRS survey analysis by jBeutler et al. (20121 
reports a value of /(0.067)(T8 (0.067) = 0.423 ± 0.055. Their 
computation is based on measuring the redshift space cor¬ 
relation function in 2D. Since the effective redshift is very 
low, their measurement is insensitive to the AP effect. 

(ii) The Sloan Digital Sky Survey (SDSS) Data Release 
7 Main Galaxy Sample analysis by Howlett et al. (20151 
reports the measurement of the two-point correlation func¬ 
tion. Using only RSD effect they obtain /(0.15)(J8(0.15) = 
0.44^0 When the AP effect is includecj^ they obtain 
/(0.15)(J8(0.15) = 0.49l°;;t 

(iii) The analysis of the Luminous Red Galaxy (LRG) 
sample in the Data Release 7 of the Sloan Digital Sky Sur¬ 
vey II by Oka et al. (20141 measures the galaxy monopole, 
quadrupole and hexadecapole power spectrum multipoles. 
They use the RSD and AP effect to report f(0.3)ag{0.3) = 
0.49 ±0.09. 

(iv) The WiggleZ analysis by Blake et al.j ( 2012| ) mea¬ 

sures the galaxy correlation function and power spectrum. 
Considering both RSD and AP effects, they report fag 
at 3 different redshifts, /(0.44)(J8(0.44) = 0.413 ± 0.080, 
/(0.6)cr8(0.6) = 0.390T0.063 and /(0.73)cr8(0.73) = 0.437± 
0.072. _ _ 

(v) The analysis of VIPERS by de la Torre et al. (20131 


® In this result the AP effect is only partially included since Qy 
and a± are set to be equal but freely vary. When both parameters 
are set free the reported value is fag = 0.63'^^' 2 ^. 
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measures the two-point correlation function. Using the RSD 
effect they report f{0.8)as{0.8) = 0.47 ± 0.08. Given the 
volume of the survey, the results are very insensitive to the 
AP effect, which is not taken into account. 


In Fig. we also plot the model predictions for 
PlancklS best-fit, Qrn = 0.308, when different theories of 
gravity are assumed. We work on the assumption of / = 
0.m{zR. The prediction for GR with a cosmological con¬ 
stant is 7 = 0.545, which is plotted in blue bands (the Icr 
limits). We also plot the predictions for two additional val¬ 
ues of 7, in red bands we plot 7 = 0.420 and in green bands 
7 = 0.680. 

Under the assumption of Planckl5+GY{,, the 7 factor de¬ 
pends on the dark energy equation of state, uj. For lo = —1 
we have the cosmological constant prediction, 7 = 0.545. 
We observe that in general all the results are in agreement 
with this prediction within la and 2cr confidence levels. For 
dark energy models with oj > —1, such as the parametrised 
Post-Friedman scalar field models, 7 > 0.545; whereas for 
to < —1, such as phantom dark energy models, 7 < 0.545. 
In the light of results of Fig. |12[ we observe that qualita¬ 
tively, redshift galaxy surveys observations slightly favour 
models with oj > —1, if ACDM -I- GR is assumed, although 
the deviation with respect w = — 1 is not very significant. 
On the other hand, if the GR condition is relaxed, 7 can 
also change. This is the case for example of the DGP model 
i Dvali, Gabadadze fc Porrati||2000 1, whose prediction for 7 
is 0.68, and it is slightly favoured with respect to GR. 


7.5 Correlated Measurements 

In this subsection we present the multivariate Gaussian 
likelihoods calculated from the DR12 data from which 
marginalised results were presented in Table We focus on 
those parameters which are of cosmological interest, such as 
fag, H{z)rs{zd) and DA{z)/rs{zd). The errors of these pa¬ 
rameters are highly correlated, as one can infer from Fig. 
Therefore, to jointly use these data, one needs to use their 
covariance matrix. 

We start by defining the data vector containing the 
cosmology parameters of interest f{z)as{z), H{z)rs{zd) (in 
lO^kms”^ units) and Da{z)/ rs{zd), 

( f{z)(^8{z) \ 

^data(^) ^ I [io®kms“^] ] . (19) 

\ DA{z)lrRzd) / 

The best-fit values of the LOWZ and CMASS samples for 
these parameters are presented in Tablej^ given a covariance 
matrix, either QPM or MD-Patchy. 

The covariance matrices of these parameters are, 

/4.1028 25.549 7.4600\ 

C'Lowf™'' = 10"® - 310.08 51.366 , (20) 

\ - - 34.912/ 

/3.7082 22.721 7.2898\ 

CgoWz = 10“" - 301.46 50.403 , (21) 

V - - 32.718/ 

/1.3424 10.597 3.7495\ 

C'cmI'sF"'' = 10"® - 179.80 34.180 , (22) 

\ - - 23.495/ 


/1.4475 11.244 4.0507\ 

C^mass = 10"' - 188.69 36.234 , (23) 

\ - - 24.698/ 

for LOWZ and CMASS samples, using the QPM or MD- 
Patchy mocks to infer the best fit parameters from the 
power spectrum multipoles measurements, as labeled. These 
covariance matrices are symmetric by construction {Cij = 
Cji) and consequently we only provide the results of half of 
the matrix. 

From these matrices, the likelihood of any cosmological 
model is given by, 

£ a exp — (D —D ) C [D —D )/2 , 

(24) 

where is the vector with the model prediction for the 

same cosmological parameters as 

In Fig. we show the ellipses which represent the like¬ 
lihood surface of Icr (Ax^ = 2.30) and 2cr (Ax^ = 6.17), 
corresponding to the covariance matrices presented above, 
for LOWZ and CMASS samples, left and right panels, re¬ 
spectively, for MD-Patchy (red lines) and for QPM mocks 
(blue lines). Each ellipse is centered in the minimum solu¬ 
tion presented in Table We observe that the differences 
between the covariances are small, and most of the shifts 
are in the minimum where they are centered, and not in the 
shape and orientation of the ellipsoid itself. Therefore, we 
conclude that the difference in covariance obtained by using 
either MD-Patchy or qpm mocks is not significant. 


8 CONCLUSIONS 


In this paper we have presented a measurement of the 
isotropic and anisotropic power spectrum relative to the 
LOS of the LOWZ and CMASS DR12 galaxy samples of 
the Baryon Oscillations Spectroscopic Survey of the Sloan 
Digital Sky Survey III. We have analysed the redshift space 
distortions in the power spectrum multipoles and the con¬ 
straints imposed on the growth factor times the amplitude 
of linear power spectrum, fag. We have also considered 
the geometrical Alcock-Paczynski effect, which allows us to 
set constrains on the angular diameter distance parameter 
DA( 2 eff)/Ts( 2 d) and the Hubble parameter p[{zes)f'sizd)- 
We study the BAO peak position on the monopole and 
quadrupole power spectrum in a companion paper which 


is released at the same time of the present paper (Gil-Marm 
et al.|2015b| ). 


In order to extract cosmological information from the 
galaxy power spectrum multipole measurements we have 


used a non-local and non-linear bias model (McDonald & 


Roy 20091 which depends on 4 parameters, bi, & 2 , b^i. and 


&3ni. Imposing that the bias is local in Lagrangian space, 5^2 
and 631,1 are set by the value of 61, and therefore only two 
free parameters are left to marginalise over, b\ and 62 - To 
model the RSD we use the TNS mod el |Taruya, Nishimichi| 
& Saito (20101; Nishimichi & Taruya (2011 1 which has been 
used in previous data releases of BOSS to describe the power 
spectrum multipoles, as well as the bispectrum monopole. 
The RSD model depends on the logarithmic rate of struc¬ 
ture growth, /, on the FoG damping parameter, ctfoG and 
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Figure 13. Likelihood surfaces of the parameters fcrg, /rs{z^) and H{z)rs{z^) extracted from the DR12 data, LOWZ sample 

(left panel) and CMASS sample (right panel), using the covariance matrices provided by MD-Patchy mocks (red ellipses) and QPM 
mocks (blue ellipses). Each ellipse correspond to Icr (Ax^ = 2.30) and 2fT (Ax^ = 6.17) confidence levels, when the 3 parameters are 
marginalised. Each type of ellipse has been centered on the minimum solution presented in Table 


on the value of the shot noise. Although / is directly related 
to the assumed given a theory of gravity, we have kept 
it free in order to test possible deviations from GR. In or¬ 
der to model the real space quantities in the TNS model, 
we have used the resumed perturbation theory at 2-loop or¬ 
der presented in Gil-Marin et al. 120121. In our model we 
have also included the geometrical AP effect, through the 
dilation parameters ay and ay, which modifies the wave 
modes parallel and perpendicular to the LOS, respectively. 
These parameters are related to the angular diameter dis¬ 
tance and the Hubble parameter, which we are also able to 
constrain. In our analysis we have fixed the shape of the lin¬ 
ear power spectrum using the fiducial cosmology but 

we have marginalised over the amplitude erg. In total, our 
galaxy redshift space power spectrum model has 8 free pa¬ 
rameters we marginalise over. 

We have tested possible systematics of our model using 
the MD-Patchy mocks simulations, which have been de¬ 
signed to reproduce up to fc ~ 0.3 hMpc“^ the power spec¬ 


trum and bispectrum of data and simulations (Kitaura et al. 


2015 companion paper). Using the mean from the 2048 sim¬ 


ulations to gain signal, we find that we are able to recover 
fag with an accuracy of < 3% for LOWZ and < 2% for 
CMASS with a fcmax ^ 0.24/iMpc“^. The systematics ob¬ 
served are significantly smaller than the statistical errors 
we measure for the DR12 BOSS data, and therefore we do 
not consider to correct the measurements, neither the cor¬ 
responding errors. 


We have computed the full covariance matrix of the 
power spectrum multipoles using two different types of 
galaxy mocks, 1000 realizations of QPM mocks and 2048 re¬ 
alizations of MD-Patchy mocks, and we have performed 
two parallel analyses using these two covariance matrices. 
Since the differences in the best fit parameters and their 
corresponding errors have been found to be small, we have 
decided to take the average among these two approaches, to 
generate a unique set of results. 

We find that for the DR12 LOWZ sample 


/(21owz)(T8(2:1owz) = 0.394 ± 0.062, DA{ziovj^)/rs{zd) = 
6.35 ± 0.19, iL(2iowz)r4zd) = (11.41 ± 0.56) 10®kms-\ 
where Ziowz = 0.32. For DR12 CMASS we find 

f(.z cmass )c^8( ■S^cniass ) — 0.444 rb 0.038, DA{Zcina3s)/l"s{zd) — 
9.42 ± 0.15, H{zcn,33s)rs{zd) = (13.92 ± 0.44) 10®kms-\ 
where Zemass = 0.57. A covariance matrix for these measure¬ 
ments was also presented. These are the main results of this 
paper and are in general agreement with previous BOSS 
DRll measurements. Furthermore, we have been able to 
reduce the previous error-bars on fas, shrinking them down 
to 15% for LOWZ sample and 8.5% for CMASS sample, 
which are the most precise measurements of this parameter 
at these redshifts at the moment, when the full AP effect 
is considered. Additionally, if we assume that the Hubble 
parameter and angular distance parameter are fixed at fidu¬ 
cial ACDM values, we find /( 2 iowz)o'g(ziowz) = 0.485±0.044 
and /(2cmass)o'8(2cmass) = 0.436 ± 0.022 for the LOWZ and 
CMASS samples, respectively. In this case the error bars 
represent a 9.1% for LOWZ and 5.0% for CMASS. 


Moreover, we have analysed the data with two addi¬ 
tional cosmological models, Qrn = 0.292 and flm ~ 0.332, 
which in this case we have changed the shape of the lin¬ 
ear power spectrum of our model, accordingly. Overall, we 
observe that both CMASS and LOWZ galaxies data is in 
agreement with the fiducial model = 0.31 consistent 
with PlancklS data: for the CMASS sample the tension is 
below la and for the LOWZ sample within 2cr. For both 
samples, when the cosmological model is changed, the value 
of fag only changes by ^ 0.5cr. 


The constraints on /(2eff)o'82eff, along with H(zeB)rd(z) 
and Da{zes)rd{zef!), will be useful in a joint analysis with 
other cosmological data sets (in particular CMB data) for 
setting stringent constraints on neutrino mass, dark energy, 
gravity, curvature as well as number of neutrino species. 
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APPENDIX A: EFFECT OF SYSTEMATIC 
WEIGHTS AT LARGE SCALES 


In this appendix we study the impact of the systematic 
weights in the power spectrum monopole and quadrupole 
of the CMASS sample. The systematic weights are de¬ 
signed to correct for fluctuations in the target density 


caused by changes in the observational efficiency (Ross 


et 


al. in prep. I. The CMASS sample presents correlations be¬ 


tween the galaxy density at large scales and the systematic 
weights are designed to correct for these variations giving 
an isotropic weighted field. However, the accuracy of these 
weights is limited and we are interested in their accuracy not 
only for the monopole, but also for the quadrupole, which 
we expect to be more sensitive to these corrections. 

In order to test this effect, we have analysed the CMASS 
data power spectrum multipoles before (pre-systematic- 
weight correction power spectrum, Pno-sys) and after the ap¬ 
plication of the systematic weights (post-systematic-weight 
correction power spectrum, Pays). We expect that these pre- 
and post-systematic-weight power spectrum multipoles con¬ 
verge to the same values at sufficiently small scales, where 
the fluctuations in the target density are not relevant. How¬ 
ever, at larger scales, where the effects of the fluctuations 
in the target density are not negligible, they will predict 


CMASS sample 



k[hMpc''] 

Figure Al. Effect of systematic weights in the CMASS data 
sample for the NGC (blue dashed lines), SGC (red dashed lines) 
and NGC-I-SGC (black solid lines). The arrows indicate the large 
scale cuts applied in order to keep those scales where the cor¬ 
rection represents less than ~ 5% of the signal. These fc-cuts 
are k = 0.02/iMpc“^ for the monopole (lower sub-panel) and 
k = 0.04/iMpc”^ for the quadrupole (upper sub-panel). 


different power spectrum amplitudes. By analyzing this dif¬ 
ference, we will estimate the percentile correction of the sys¬ 
tematic weights on the power spectrum multipoles. Setting 
a limit of ~ 5% correction, we define a large scale cutoff for 
the power spectrum multipoles. 

This is displayed in Fig. ig where the ratio between the 
data power spectrum multipoles is shown for the monopole 
(lower sub-panel) and quadrupole (upper sub-panel). The 
different colour lines show this effect for the NGC, SGC, 
and NGC-bSGC as labeled. 

As expected, the effect of the systematic weights is rel¬ 
evant only at large scales, where they suppress spurious cor¬ 
relations and has a higher impact on the quadrupole respect 
to the monopole. Setting an accuracy limit of ~ 5% correc¬ 
tion, we discard those scales where the systematic weight 
correction on the power spectrum multipoles exceeds this 
limit. Furthermore, we can compute the for the differ¬ 
ence between the weighted and unweighted and P^^\ 
This is x^ = DC~^D^, where D = Psya — Pno-ays. We ob¬ 
tain that for fcmin = 0.02feMpc“^ for the monopole and 
fcmin = 0.04/iMpc“^ for the quadrupole, (as indicated by 
the black arrows) and femax = 0.25/iMpc”^ for both statis¬ 
tic, the reduced x^ is 0.044. This means that by applying 
the weights we are maximally correcting by just ~ 0.2a 
(0.2 ~ VO.044) in some measured parameter. This set a 
maximal impact of the known systematic without the need 
of defining any particular model. 

Therefore, in this paper we will only consider for our 
analyses those scales smaller than k = 0.02 /iMpc”^ for the 
monopole and k = 0.04ftMpc“^ in the quadrupole, as indi¬ 
cated by the black arrows. 
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DR12 dataset QPM MD-Patchy 

LOWZ 0.0157 0.0418 0.0148 

CMASS 0.0528 0.0529 0.0349 

Table Bl. Fraction of collided galaxies, fcg = — 

1]/ u)f(.(xi), for the data, for the QPM and MD-Patchy mocks, 

for the CMASS and LOWZ samples. 


APPENDIX B: EFFECT OF FIBER 
COLLISIONS WEIGHTS ON THE BEST-FIT 
PARAMETER ESTIMATION 

In this appendix we study the effect of the fiber collision 
weights on the power spectrum multipoles, and more pre¬ 
cisely on the fas measurement according to the model pre¬ 
sented in As described in s the fiber collision weights 
are included in order to account for those galaxy pairs that 
are too close to each other (< 62") to put two fibre detec¬ 
tors. The fraction of collided galaxies, fcg = — 

1)]/ 5])^ Wfc(xi), for the data is presented in the first column 
of Table |BTj for the CMASS and LOWZ samples. Since the 
number density of galaxies is higher in the CMASS sample, 
the fraction of collided galaxies is also higher in this sam¬ 
ple. The second and the third column of Table |BT] present 
the fraction of collided galaxies found in the QPM and MD- 
Patchy mocks, respectively. For the CMASS sample, the 
MD-Patchy mocks present a smaller value of fcg respect to 
the data and QPM mocks, which both are in close agreement. 
This is due to a limitation in the resolution of substructure 
inside MD-Patchy haloes. This limitation will be solved in 
future versions of the MD-Patchy mocks. In the LOWZ 
sample, both data and MD-Patchy mocks agree well with 
the value of fcg, because for the number density of galax¬ 
ies of this sample, the resolution of substructure was not a 
limiting factor. However, we see that QPM mocks present a 
higher value of fcg with respect to the data. This is due to 
the version of QPM mocks used in this paper matched an 
old LOWZ catalogue, where galaxies with previously known 
redshifts were subsample to match the BOSS close-pairs se¬ 
lection. This subsample has now been discontinued. As for 
the MD-Patchy mocks, this will be fixed in future releases. 

In order to test the effect of the fiber collisions in the 
parameter estimation, we focus on the QPM mocks for the 
CMASS sample, which has the higher value of fcg amongst 
all the cases. We measure the monopole and quadrupole of 
1000 mock realizations and take their average in order to 
gain signal to noise. We consider the two following selection 
of galaxies 

(i) We treat the galaxies as in a real survey. When two 
or more galaxies present an angular separation of ^ 62", we 
weight one of them by the number of galaxies within the 
^ 62" angular radius, and remove the others. This mimics 
what it is done with the real dataset 

(ii) We consider all the galaxies resolved in the mocks 
and weight them equally. This is the ideal case we would 
have if all the targeted galaxies in the survey were analysed 
spectroscopically. 

The case (ii) has the correct clustering and the correct 
anisotropic signal. Comparing both will show how the colli¬ 


sion weights affect the power spectrum multipoles and the 
estimation of parameters. 

When we compare the power spectrum monopole (and 
therefore the isotropic clustering) of cases (i) and (ii), we see 
that the effects of the collision weights are degenerate with 
the amplitude of the shot noise parameter, Anoiae, which 
we treat as a free nuisance parameter. This was already re¬ 


ported in Gil-Marin et al. 12015a I for the DRll sample using 
PThalos mocks ( [Manera et al. 20131. The effect of the fiber 
collision weights in the power spectrum quadrupole (i.e. in 
the anisotropic clustering) is more complex. We know that 
a fraction of the angular close pairs correspond to galaxies 
which share the same dark matter host halo. Using the re¬ 
gions where we have a superposition of plates, we have esti¬ 
mated that this fraction is about ~ 60%. On the other hand, 
~ 40% of the angular close pairs correspond to galaxies that 
happens to share a similar LOS, but that are actually rea¬ 
sonably uncorrelated. By applying the fiber collision weights 
to galaxies that are actual close pairs, we are removing signal 
in the direction of the LOS with respect to the signal in the 
transverse direction, which is not modified. Therefore, by ap¬ 
plying the collision weights, the anisotropic power spectrum 
is affected, which potentially can alter the estimation of fas- 
This is shown in the left panel of Fig. |B1| where the black 
solid line show the fractional change in the power spectrum 
quadrupole of case (i) respect to case (ii), for QPM mocks. 
We see that this change is sub-percent at large scales, but 
rapidly grows as we go to smaller scales, reaching ~ 10% at 
k ~ 0.2 /iMpc“^ and ~ 20% at A: ~ 0.25 hMpc“^. In order to 
test the impact of this change in the fitted parameters, we fit 
the RSD bias model presented in [J^to the power spectrum 
monopole and quadrupole of cases (i) and (ii) and compare 
them, using the scales 0.02hMpc“^ ^ A: ^ kmiix for the 
monopole and 0.04/iMpc“'^ ^ A: ^ Aimax for the quadrupole, 
as we did for the data, varying A:niax from 0 .15 / iMpc~^ to 
0.24ZiMpc”'^. The blue dashed line of Fig. ^ shows the 
ratio of the best-fit models to power spectrum multipoles 
of cases (i) and (ii). We see that the difference among the 
models describes well the differences in the power spectrum 
quadrupole produced by the fiber collision weights. In the 
right panel of Fig. |B1| the ratio between the best-fit pa¬ 
rameters, {fel(T8,fe2Crg,/(J8,CTPoG, Anoise,a||,a:Y} of theSB twO 
models are shown. Both top and bottom panels display the 
same ratio of parameters, but at different ranges, for clarity. 
We observe that the change in the monopole and quadrupole 
due to the fiber collision weights is absorbed mainly by Anoise 
and b 2 as, which changes of order of 20% — 100% in the fiber 
collision case respect to case (i). The FoG damping param¬ 
eter, cTpoG changes by about 2.5% and does not present any 
significant change with the minimum scale of the fit. The AP 
parameters, ay and a± do not present any significant change 
due to the effect of fiber collisions. Finally fas and bias are 
modified by about ~ 0.5%, and show no dependence with 
the minimum scale of the fit. In Table [B2l we summarise 
these results and compare them with the statistical errors 
for the CMASS sample using Amax = 0.24/iMpc“^. 


We conclude that the effect of fiber collision are ab¬ 
sorbed chiefly by Anoise and b 2 as- The impact of the fiber 
collision on fas and bias is a sub-percent. In this case, fas 
tend to be overestimated with respect to the ideal case where 
all the galaxies were considered, whereas bias tend to be un- 
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Figure Bl. Effect of fiber collision weights in the CMASS sample, where the fraction of collided galaxies is 0.0529. The left panel display 
the ratio between power spectrum quadrupole of the cases (i) and (ii) (see text), where the fiber collisions are applied, Collision Weights? 
and where all the galaxies are considered, AiiGaiaxies- The black solid line is the ratio of the measured quadrupole of 1000 realizations 
of the QPM mocks, and the blue dashed line the ratio between the best-fit models to the cases (i) and (ii) at fcmax = 0.24/iMpc“^. The 
right panel displays the ratio between the best-fit parameters of the best-fit models of the cases (i) and (ii). Top and bottom panel display 
the ratio at different scales. 



QPM mocks (CMASS sample) 



kmax [hMpc''] 


Parameter 

statistical error [%] 

systematic error due to fc [%] 

/<78 

8.5% 

0.5% 

"II 

3.1% 

0 .1% 

a_L 

1 .6% 

< 0.1% 

bias 

1.7% 

0.5% 

fees 

130% 

70% 

^FoG 

9.5% 

1.5% 

A 

-^noise 

175% 

10 % 


Table B2. Statistical and systematic errors caused by the effect 
of fiber collisions in the free parameters of the model, for the 
CMASS sample at kmax = 0.24 hMpc“^. For the cosmological 
parameters, /erg, q;|| and a± the systematic errors caused by the 
fiber collisions are much smaller than the statistical errors of the 
data for the CMASS sample. 


derestimated. Finally the AP parameters present changes of 
order ~ 0.1% for ay and 0.1% for a±. 

Since the Icr statistical errors of these parameters are 
about 10 times larger than this systematic shift, we do not 
correct our result by this effect. Because the fraction of col¬ 
lided galaxies in the CMASS sample is higher than in the 
LOWZ sample, we expect that these changes are also negli¬ 
gible in the LOWZ sample. 
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